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1 | Inleiding 


De afgelopen maanden zijn weer behoorlijk druk geweest voor het CAN. De toewijzing van 
het Research Institute for Applications of Computer Algebra (RIACA), waarvan in de vorige 
Nieuwsbrief al sprake was, bleek toch vollediger uit te vallen dan we op dat moment dachten. 

Dat betekende dat in korte tijd een hoop plannen gemaakt moesten worden, zowel op 
wetenschappelijk als op organisatorisch terrein. Hoewel het wachten op dit moment nog is op 
de uiteindelijke toekenning van de gelden voor 1993 en er nog geen concrete stappen genomen 
kunnen worden, betekent de uitbreiding van de CAN-activiteiten onder meer dat de huidige 
locatie bij het CWI verlaten moeten worden: er is geen ruimte meer voor uitbreiding. Het 
laat zich aanzien dat in dat geval de huisvesting plaats zal vinden in het Matrix gebouw, op 
het WCW terrein, dat onlangs officiéel geopend is. 

Daarnaast heeft CAN als codrdinator een HCMP-aanvraag voorbereid en ingediend voor 
een Europees netwerk voor symbolisch en algebraisch rekenen. Doel is te komen tot een 
Europees platform voor computeralgebra vanwaaruit diverse activiteiten gecoordineerd en 
ondersteund worden. Onderzoeksgroepen uit acht landen nemen deel aan dit project. 

De gewone CAN activiteiten gaan ondertussen gewoon door, na de personeelsuitbreiding 
zelfs met verdubbelde kracht. Veel van het werk dat de afgelopen tijd verricht is, was gericht 
op de interne organisatie en niet direct zichtbaar, maar we hopen uiteraard dat de gevolgen 
zichtbaar zullen worden door de manier van werken binnen het Expertisecentrum. 


Jan Sanders 
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2| News from the CAN Expertise Centre 


2.1 New Computing Facilities for Computer Algebra 


Recently, new central servers at the CAN Expertise Centre have been configured for computer 
algebra systems and have been made available via computer networks. The new equipment 
replaces the previous CAN computers CANA and CANB. It consists of 


e two Silicon Graphics Indigo Servers, called CANC and CAND, with a R4000 processor, 
64 MB RAM and 1 GB disk space. 


e a multi-processing Silicon Graphics Challenge L server with two R4400 processors, 128 
MB RAM and 2 GB disk space. This machine is called CANE. 


Typical computer algebra tests with Maple, REDUCE and FORM show that the new 
equipment is five to eight times as fast as the old SUN4/280-computers. 

The usage of the central machines is currently set up by institutional contracts and indi- 
vidual accounts. A computer algebra service arrangement with institutions that have many 
potential computer algebra users offers them not only the access to the CAN computers, but 
also an extensive support programme from the Expertise Centre, like programming help and 
sessions for users. 

Beside the computer algebra service arrangement, the possibility of a strictly-personal 
login is still possible. A registration form is included in this newsletter. 


2.2 Upcoming CAN Courses and Events 


The Computer Algebra Netherlands Expertise Centre organises in this Fall several courses in 
the use of computer algebra systems. Scheduled are introductory courses and programming 
courses in Axiom, Maple and Mathematica. The dates are: 








Date Course 
September 30 | Introduction to Maple 
October 1 Programming in Maple 


November 1 Introduction to Mathematica 
November 2. | Programming in Mathematica 
December 15 | Introduction to Axiom 














Furthermore, a computer algebra symposium on the occasion of the opening of the Re- 
search Institute in Applications of Computer Algebra (RIACA) is scheduled in October. 


2.3. Distribution of Mathematical Software 


The product range of commercial mathematical software which is available through CAN has 
recently been extended to LOCBIF, MACSYMA, and MATHTENSOR. Elsewhere in this 
newsletter you can find specific product information. 
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2.4 Let Me Introduce Myself 


In keeping up the nice tradition of this newsletter I would like to introduce myself as the 
new system administrator for CAN computers starting May 1993. My name is Alena Lacova 
and I was born in Slovakia (previously part of Czechoslovakia). I have been interested in 
mathematics and physics for a long time. I enjoy solving puzzles and took part in various 
competitions. In 1982 I enrolled in the Faculty of Mathematics and Physics of the Comenius 
University in Bratislava, Slovakia and received a M.Sc. in nuclear physics in 1987. Afterwards 
I accepted a study position at the Department of Nuclear Physics and Technology of Slovak 
Technical University in Bratislava where I was responsible for developing educational software 
and installing software for the experiments I was involved in. In 1989 I joined the computer 
group at National Institute for Nuclear and High Energy Physics in Amsterdam, as a UNIX 
system administrator and programmer. I was also providing support to users with PostScript 
problems. Yet recently I found an occupation which is more challenging and satisfying better 
my interests. It is the combination of my interests in computers and mathematics, i.e. the 
system administrator maintaining the CAN and RIACA computing facilities and giving the 
user support in UNIX, PostScript and partially also CA related problems. So if you have any 
problems when using CAN computers or with PostScript or you need a new account, I will 
be pleased to help you. You can send e-mail to alena@can.nl or phone me on 020-5926050. 


Alena Lacova 
alena@can.nl 


In order to cope with the strong growth of CAN more personnel was needed to manpower 
the increase. This success shows that the promotion of CA systems has been well taken care 
of. I have been contracted to ensure progress will continue and that the need for governmental 
support will decrease. My task will be to market CAN and to run the ’shop’, that is making 
the software available to you. 

My name is Dick Verkerk. I was born in Losser, which lies in the east of the Netherlands. 
A great part of my life I lived in Germany. After my study International Marketing in 
Groningen I have worked for various organizations on short term contracts. Together with 
Els Willems I will be your contact for the availability of software. Technical questions are the 
expertise of André Heck and Leendert van Gastel. You can send e-mail to verkerk@can.nl or 
phone me on 020-5926051. 


Dick Verkerk 
verkerk@can.nl 


Leendert van Gastel, André Heck 
CAN Expertise Centre 
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3 | From the Annual Report 1992 


3.1 General 


In mathematics symbolic computation on a computer, called “computer algebra” consider- 
ably increases the computing power at hand, and can strongly improve the working efficiency. 
Series of symbolic instructions can be formulated in a programme, and thus routine manip- 
ulations can be handed over to the computer. 

After the implementation of numeric computation and graphics, symbolic computation 
was the missing link in the chain of parts of mathematics that could be automated. With the 
development of the algorithms and the growth in power and speed of the computer, computer 
algebra became available to a large audience. 

It is now possible to work with “general purpose systems” that offer symbolic, numeric 
and graphical features, and that support a programming language, which is close to the 
original language in which the problems are posed. These new possibilities lead to a different 
way of viewing mathematics and computers. 

Also a different point of view is needed on the education of mathematics. Software systems 
that have enough built-in knowledge to pass the secondary school examinations, have a large 
impact on the didactics of the sciences. This impact is comparable to that of the pocket 
calculator. 


3.1.1 Goal 


The goal of the CAN Foundation is the stimulation and support of computer algebra in 
research, education and development. With computer algebra we mean the design, analysis, 
implementation and usage of algebraic algorithms in the widest sense. Among the means to 
arrive at this goal are: 


e the development of knowledge and expertise in the field of computer algebra; 
e a broad range of computing facilities for the use of computer algebra; 


e the coordination of activities like courses and seminars on this field. 


In 1989 an Expertise Centre was established to embody all the activities. It is one 
in the row of Expertise Centres that have as their task to provide disciplinary knowledge 
with use of advanced information and communication technology. There are about twenty of 
these Expertise Centres that are recognised by NWO (Netherlands Organisation for Scientific 
Research). Most of them have received an establishment grant. 


3.1.2 Activities in 1992 


In a nutshell we present the activities in 1992 of CAN. Elsewhere in this annual report, these 
activities will be dealt with in greater detail. 


10 Nieuwsbrief 


1. Computing facilities have been made available for the special purpose of computer 
algebra: workstations located at various research groups that work intensively with 
computer algebra, and central servers that are configured for symbolic computation. 
With these computing facilities all packages that are of interest to mathematicians, 
physicists, engineers and others, can be reached via Internet or SURF net. 


2. Information and support for the usage of computer algebra. Here the usage of infor- 
mation technology tools is prominent. Electronic software libraries and an electronic 
FTP archive are maintained. The help desk uses e-mail as much as possible, though 
telephone and fax are used as well. A freely accessible information system for all kinds 
of data that concern computer algebra, is in development. Only the newsletter that 
comes out every half year, is printed on paper. In this case, paper remains still the best 
medium. 


3. Coordination and stimulation of the use of computer algebra in education via a working 
group, named CONCAVO. The aim is to let education profit from the new developments 
by distributing professional and easy-to-use mathematical software, and by furthering 
cooperation by exchanging material and experiences. There is a preprint series to 
facilitate the distribution of material. 


4. Distribution of software systems developed by CAN or others. CAN offers two special 
purpose systems LiE and FORM to an international market. Moreover, general purpose 
systems like DERIVE, Maple and Mathematica are in the product range. 


5. Organisation of seminars and courses. The susceptibility for a symbolic computational 
approach of researchers from industry and business can be extended, for example by 
giving supporting the mathematical formulation of their problems. This was the aim 
of the SCAFI project (Studies of Computer Algebra For Industry). In 1992 a second 
seminar for SCAFI was organised by CAN. Furthermore an introductory course in 
Mathematica was held in spring. The overwhelming response shows that there is a 
large interest in this way of supplying expertise. 


The activities of CAN are directed to three main fields: 


e Scientific research, with an emphasis on the sciences; 
e Education, ranging from secondary school to university; 


e Business: industrial research, engineers and finance. 


3.1.3. Trends and Perspectives 


The interest in computer algebra and the usage of these systems has grown substantially in 
1992. Already on a local level the knowledge about the general purpose systems starts to 
grow. Furthermore, in the last few years, the aS power available for research groups 
has grown tremendously. 

For the activities of CAN, these developments brought a shift: from hardware towards 
software, from computing power towards information and support. The initiative to provide 
research groups with workstations to make automated symbolic computation possible reached 
its goal prior to expectation, because of the growth in available computers and software. The 
central computing facilities, however, still are of great value to the scientific research in the 
Netherlands, because of their configuration as “symbol crunchers” and the special software 
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that they carry. CAN remains unique as a centre where all kinds of computer algebra software 
can be found. 


The growing interest makes computer algebra a fertile soil. Quite often activities initi- 
ated by CAN have been taken up by others. Publishers showed interest in the CONCAVO 
preprints, two of them are going to be published. Moreover several publishers showed interest 
in the CAN newsletter. 

The seminar SCAFI attracted such a lot of international attention that apart from the 
seminar in Amsterdam, both in Bath (UK) and in Brussels seminars with the name SCAFI 
were held. 


International cooperation is of utmost importance for a national centre if it is to function 
in an optimal way. 


e CAN has put a lot of effort in making also the SCAFI seminars in Bath and Brussels 
a success. 


e CAN has given assistance to Wolfram Research in organising the European Mathematica 
Conference in Rotterdam. This attracted more than 300 visitors. 


e An initiative has been launched to set up an electronic information system together 
with partners in the United States, Denmark and Germany. 


e CAN has functioned as coordinator for a proposal for a European network for Symbolic 
and Algebraic Computation within the framework of the EEC Human Capital and 
Mobility programme. 


Right at the end of 1992 it was announced that the proposal for RIACA was accepted for 
government funding. RIACA, a collaborative initiative of SMC (Dutch Mathematical Centre 
Foundation), RISC-Linz (the well-known Research Institute for Symbolic Computation) from 
Austria, and CAN. It will comprise a research institute of a modest size, but with a staff of 
international top level. On project base applications of computer algebra will be worked out, 
with attention to the theoretical background and to the technical implementation. RIACA 
can serve as a bridge between theory of mathematics and computer science, and every day 
practice of computer and software. 


Though the income of CAN clearly shows growth, the finances of the Expertise Centre 
are not without its difficulties. Also the support of SMC does not completely fill the gaps. 
Finances also worried other NWO Expertise Centres, and the Dutch Minister of Education 
and Science has created some possibilities for funding in cooperation with NWO. CAN hopes 
that with this support it will be able to match expenses and income from 1995 onwards. 


3.2 Addendum 


Please, contact the CAN Expertise Centre if you are interested in the complete Annual Report 
1992. 
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4| CONCAVO 


CONCAVO is een werkgroep ingesteld door CAN voor de uitwisseling van ervaringen met 
het gebruik van computeralgebra in het onderwijs. Dientengevolge zijn er drie niet strikt 
gescheiden werkterreinen: het universitair onderwijs, het hoger beroeps onderwijs en het 
vwo. Op ieder van deze terreinen zijn er, met wisselende intensiviteit, activiteiten. Het 
volgende beoogt daar een klein overzicht van te geven. 


4.1 Het universitair onderwijs 


Bij de start van CONCAVO waren er wat contacten tussen personen die actief bezig waren 
computeralgebra in het onderwijs te betrekken. Een gespreksgroep binnen CONCAVO is 
echter niet van de grond gekomen; incidentele persoonlijke contacten bestaan nog wel. 

Bij het universitair onderwijs kan computeralgebra als zodanig een doel zijn of het kan 
als stuk gereedschap of anderszins een plaats hebben in andere vakken. 

Onlangs is een kleine telefonische inventarisatie gedaan naar het gebruik ervan in het 
onderwijs. Bij vrijwel alle universiteiten blijkt daar belangstelling voor te bestaan. Bij 
alle universiteiten is computeralgebra beschikbaar. Hier en daar worden er bij bestaande 
vakken praktika gegeven. Ook wordt op sommige plaatsen onderwijs in het vak computeral- 
gebra of in bestaande pakketten, zoals MAPLE of MATHEMATICA, gegeven. Vaak hangen 
de activiteiten sterk af van de aanwezigheid van geinteresseerde stafleden. In het kader van 
CONCAVO zal geprobeerd worden deze geinteresseerden met elkaar in contact te brengen. In 
eerste instantie zal van iedere universiteit een (informele) vertegenwoordiger gezocht worden; 
op dit moment moet nog bij twee universiteiten iemand gevonden worden. Deze vertegen- 
woordigers zullen binnenkort uitgenodigd worden voor een gesprek op een centrale plaats in 
Nederland, waar dan afspraken gemaakt kunnen worden voor eventuele verdere activiteiten. 


4.2 Het hoger beroeps- en technisch onderwijs 


Er is bij het hoger beroeps onderwijs veel aandacht voor computeralgebra. De gesprekken 
binnen CONCAVO vinden op dit moment echter alleen plaats tussen diegenen die er twee jaar 
geleden ook bij betrokken waren. Op 16 januari was er een bijeenkomst in Tilburg waar o.a. de 
situatie rond CONCAVO en het hoger beroeps- en technisch onderwijs besproken is. Het lijkt 
zinvol ook anderen bij de CONCAVO-gesprekken te gaan betrekken. Een goede gelegenheid 
daarvoor is de computeralgebra dag voor het HBO, die op 12 mei gemeenschappelijk door 
CAN en de Noordelijke Hogeschool Leeuwarden georganiseerd gaat worden. 


4.3 Het vwo 


De eerste activiteiten van CONCAVO op dit terrein lagen op het produceren van werkbladen 
voor het vwo met het pakket DERIVE. Deze activiteit werd uitgevoerd door Paul Drijvers, 
die het resultaat inmiddels heeft gepubliceerd. Bij de deelnemers is belangstelling om verder 
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te gaan. In de onderwijswereld is nu veel belangstelling voor de grafische calculator, waarmee 
veel eenvoudiger dan met de hand functieonderzoek gedaan kan worden. Het effect van 
computeralgebra zal nog veel groter zijn. Zodra een pakket als DERIVE op een goedkope 
pocketcomputer beschikbaar zal zijn (en dat is slechts een kwestie van tijd), kunnen de 
leerlingen efficienter formulemanipulaties met deze pocketcomputer doen dan met de hand. 
De vraag is hoe het onderwijs op deze ontwikkelingen in zal gaan spelen. Er is hierover contact 
geweest met de Nederlandse Vereniging van Wiskunde Leraren. Op een drietal regionale 
bijeenkomsten, in Amsterdam, Zwolle en Eindhoven is aan dit onderwerp aandacht besteed 
en dat heeft geleid tot een aantal nieuwe deelnemers voor de werkgroep over computeralgebra 
in het vwo. Binnenkort zal‘een bijeenkomst belegd worden om de toekomstige activiteiten te 
gaan bespreken. 

Overigens dienen bij de activiteiten van CONCAVO op dit terrein ook de lerarenoplei- 
dingen betrokken te worden. 

Nadere informatie over de activiteiten van CONCAVO zijn te verkrijgen bij CAN, en bij 


Fred Simons 

Faculteit Wiskunde en Informatica 
Technische Universiteit Eindhoven 
Postbus 513 : 
5600 MB Eindhoven 

tel. 040 - 474400 
wsgbfs@win.tue.nl 


Ernic Kamerich 

Faculteit Wiskunde en Informatica 
Katholieke Universiteit Nijmegen 
Toernooiveld 
6525 ED Nijmegen 

tel. 080 - 653220 / 652986 
ernic@sci.kun.nl 
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5 | New Mathematical Software 


5.1 LOCBIF 


In this short note you’ll learn about the general features of LOCBIF, the authors of LOCBIF, 
system requirements for using LOCBIF, and pricing. 


5.1.1 What is LOCBIF ? 


LOCBIF is a new powerful tool for bifurcation analysis of ordinary differential equations 
(ODEs) which depend upon parameters. It allows you to explore interactively the existence 
and stability of equilibria in dynamical models! . 

LOCBIF is based on a continuation procedure for relevant local bifurcation curves up 
to codimension three. Projection of these curves onto the parameter space determines the 
boundaries of the equilibrium existence and stability. 

LOCBIF combines modern results on normal forms and local bifurcations with interac- 
tive computer graphics, resulting in a unique integrated environment for ODE model analysis. 

LOCBIF plots two dimensional projections of the bifurcation curves during the compu- 
tation which may be stepwise or automatic. 

LOCBIF maintains an archive of ODEs and allows you to specify a new system of ODEs 
during a session in a simple Pascal-like language. It compiles the right hand sides of ODEs 
by means of an on-line compiler. Computed curves may be stored on a disk. Stored curves 
may be plotted and used in further computations. 

LOCBIF allows you to study dynamical systems with up to ten phase variables and ten 
parameters. 

LOCBIF is designed for students, teachers and scientists. 


5.1.2 Authors of LOCBIF 


LOCBIF combines numerical algorithms developed by A. Khibnik and E. Nikolaev with 
modern interface designed by Yu. Kuznetsov and V. Levitin. 

The numerics of LOCBIF are based on an improved version of the LINLBF code devel- 
oped by A. Khibnik. The new continuation code BEETLE used in LOCBIF was developed 
by E. Nikolaev and A. Khibnik. 

The general dialog scenario and graphics window interface of LOCBIF were designed 
and implemented by Yu. Kuznetsov. 

The built-in RHS editor and on-line compiler, as well as a manager for system archive, 
were designed and written by V. Levitin. 

As in any large project, LOCBIF is more than simple sum of the parts mentioned above. 
A number of useful features of LOCBIF are the result of close collaboration of all the authors. 





‘Three additional versions of LOCBIF are also available: for analysis of fixed points and periodic or- 
bits of iterated maps, periodic solutions of periodically forced ODEs and periodic solutions (limit cycles) of 
autonomous ODEs. These versions are presented in Appendix C. 
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For general reference concerning LOCBIF see (Khibnik, Kuznetsov, Levitin, Nikolaev, 
1992). An earlier presentation of the numerical algorithms implemented in LOCBIF can be 
found in (Khibnik, 1990a,b). 


The permanent address of the authors of the LOCBIF: 


Institute of Mathematical Problems of Biology, 
Russian Academy of Sciences, 

Pushchino, Moscow Region, 

142292 Russia. 


E-mail: impb@venus.iteb.serpukhov.su 


5.1.3 System requirements 


LOCBIF requires an IBM PC/XT, AT, PS/2 or compatible microcomputer running DOS 
version 3.0 or higher. LOCBIF automatically supports monochrome and color EGA, VGA 
and Hercules monochrome graphics display adapters. It requires 640K of RAM with no less 
than 500K free. A math coprocessor is not required but LOCBIF will use it if one is present. 
The math coprocessor will greatly decrease computation time and is strongly recommended. 

LOCBIF can be run from a double or high density floppy disk but it will work more 
efficiently if all its files are on a hard disk. 

If your computer has more than 120K of extended memory, LOCBIF can use a virtual 
(RAM) drive installed in this memory to store working files. 


5.1.4 Pricing 
LOCBIF is a commercial product and is available from the CAN Expertise Centre. 


License Form Price 
Single User Dfi. 500 
Multi User license (allows use at several machines at one site) Df. 1,500 


Campus License (as Multi User and allows home-use by employees) Df. 2,500 


5.2 PC Macsyma 


In the Fall of 1992 Macsyma Inc. began shipping PC Macsyma 417.125 for 386 and 486 PCs. 
This product is based on Macsyma 417.100 which has been available on several UNIX plat- 
forms for over a year. PC Macsyma 417.125 requires Microsoft Windows 3.1 and runs under 
MS-DOS 5.0. Since April 1993 Macsyma is available from the CAN Expertise Centre as well. 


Information about the product is available from the CAN Expertise Centre by phone, mail, 
fax, or email. This info includes a “Highlights of Macsyma” document, a “System Require- 
ments” document, an order form, and a listing of the major enhancements made to Macsyma 
since its licensing from M.LT. 


A free demo diskette for the product is also available demonstrating its math, graph- 
ics, and help facilities. 


New Features 
Besides some bug fixes to Macsyma 417.100, some of the important new features of PC 
Macsyma 417.125 are: 


o It is a good citizen of Windows 3.1, working in enhanced mode. 
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o It has much better performance and/or capability than the current 415.25 PC release in 
many areas including integer factoring and certain multivariate polynomial factorizations, 
prime number testing, trig special angle simplification, double precision floating point 
computations, ode solving (see below), etc. 


o It has much enhanced capabilities for 2D and 3D graphics, including parametric plotting 
of surfaces, controls for rotate, zoom, colors, etc. There is a viewer enabling looking inside 
surfaces. 


o It has a new input editor. 
o It has a new hypertext Help facility. 


o It has a new “fancy” mode for the display and hardcopy of results (Macsyma D-lines) 
featuring variable width fonts (user’s choice as to which and what size), smaller sized 
exponents and subscripts, Greek characters, appropriately sized parentheses, integral, 
summation, and product signs, etc. 


o It has yet another package (named ODELIN2) for solving ODEs, designed and imple- 
mented by Prof. Shunro Watanabe of the Department of Mathematics, Tsuda College, 
Tokyo, Japan. The ODELIN2 package solves 2nd order LINear ODEs using a technique 
based on the hypergeometric equation. Prof. Watanabe says that the program can solve 
approx. 96% of the linear ODEs of 2nd order out of almost 500 applicable equations in 
Kamke’s book. 


o Its integration capability has been improved in a few areas over Macsyma 417.100. For 
example: 


* improved handling for some trig integrals e.g. integrate (x/sin(x)?, x); > log(sin(x)) — 


x/tan(z) 
* handling some integrands involving abs e.g. integrate(xabs(ar + b),z);— (ax + 
b) (2ax — b) abs(az + b)/ (6a) 


* generalized zeta integrals e.g. integrate (x°/(exp(zx) — 1), z, 0, inf); > 14/15 


o It includes a new basic linear algebra package providing LDU decomposition, LU decom- 
position, Cholesky decomposition, support for Toeplitz and Hankel matrices, etc. 


System Configuration 

The basic hardware: a 80386 or 80486-based PC, IBM compatible. 

Currently, a math coprocessor is required. (A 80387 is required for use with a 80386 CPU. A 
80487 is required for use with a 80486SX CPU. The 80486DX includes a math coprocessor.) 
RAM: 8 megabytes minimum, 12 megabytes is recommended. 

Hard disk: 20 megabyte minimum free space must be available. 

Graphics adapter: Super VGA, VGA, EGA (or any adapter supported by MS-Windows 3.1). 
Printers: Compatible with any printer supported by MS-Windows, which includes dot ma- 
trix, Postscript, HP Laser Jet, HP Desk Jet. MS-Paint and MS-Write can make hardcopy 
bitmaps if desired. 

Network software: Works with most PC network software. A floating license manager is 
available which works with Novell networks. 

Currently, the product ships on 3.5” diskettes. 


For more complete specifications, request the "System Requirements” document. 
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5.3 MathTensor 


5.3.1 Hardware Requirements 


MathTensor requires approximately 2 megabytes of disk space. It is recommended that a 
workstation have at least 8 megabytes of RAM memory. The Macintosh version runs on 
machines with 8 megabytes or more. DOS versions require at least 5 megabytes of memory, 
but more is highly recommended. The Windows 3.x version works with 8 megabyes of RAM 
or more. 


5.3.2 Software Requirements 


MathTensor requires Mathematica 1.2, 2.0, or greater. 


5.3.3 Version 


The current version of MathTensor is 2.1.5. 


5.3.4 Contact 
MathTensor is available from 


MathSolutions, Inc. 
P.O. Box 16175, Chapel Hill, NC 27516 
Telephone/ Answering Machine/FAX: 919-967-9853, Email: mathtensor@wri.com 


and 


CAN Expertise Centre 

P.O. Box 4079, 1009 AB Amsterdam 

The Netherlands 

Telephone: +31 20 5926050, Fax: +31 20 5924199, Email: can@can.nl 


5.3.5 Authors 
Leonard Parker, Ph.D. and Steven M. Christensen, Ph.D., founders of MathSolutions, Inc. 


5.3.6 General Information 


Tensor analysis is extensively used in applications in physics, mathematics, engineering, and 
many other areas of scientific research. Problems involving tensors often are extraordinarily 
large and can be some of the most difficult computations in all of science. Equations with 
thousands of terms are common and can only be manipulated by computer mathematics 
systems like Mathematica. 

MathTensor is the largest Mathematica package yet developed outside of Wolfram Re- 
search. It adds over 250 new functions and objects to Mathematica to give the user both 
elementary and advanced tensor-analysis functionality. MathTensor is a general tool for 
handling both indicial and concrete tensor indices. Standard objects like Riemann tensor, 
Ricci tensor, metric and others are built into the system along with common functions like 
the covariant derivative, index commutation, raising and lowering of indices, and various 
differential forms operations. 

MathTensor provides commands for simplifying and manipulating tensor expressions, as 
well as a knowledge base of transformation rules and definitions required for dealing with some 
of the more important tensors. MathTensor is designed to work along with the functions of 


CAN 19 


Mathematica to provide users with the functions and objects they need to devise their own 
custom tensor analysis programs. MathTensor provides most of the basic structures needed 
for doing tensor computations and for programming new functions. As MathTensor is used 
and special applications are developed, they will be added to future versions of MathTensor. 

MathTensor has been under development by Leonard Parker and Steven M. Christensen 
since the first alpha test release of Mathematica. It contains over 25,000 lines of Mathemat- 
ica code contained in nearly 100 files totalling approximately 2.0 Megabytes of disk space. 
MathTensor will run on any machine that runs Mathematica and has sufficient RAM memory 
(generally 8 Megabytes or more) and disk space for file storage and swap. MathTensor runs 
under versions 1.2 and 2.X of Mathematica. 


5.3.7 MathTensor Examples 


In the Computer Algebra ftp-archive on canb.can.nl are two information files that contain 
examples located in ftp/pub/mathematica/MathTensor. 


5.4 New Features in Mathematica Version 2.2 


Wolfram Research Inc. and the CAN Expertise Centre are proud to announce the availability 
of Mathematica Version 2.2. Users with a Math Plus contract will receive the upgrade as 
soon as possible. In this section we give a short list of the most significant enhancements in 
Mathematica Version 2.2. A long listing with examples can be obtained by anonymous ftp 
from canb.can.nl throug the file ftp/pub/mathematica/Misc/new in _mathematica_2.2. 


5.4.1 Computational Enhancements 


e The numerical solution of sparse linear systems of equations has been made much faster. 


Symbolic definite integration now includes tests for nonintegrable singularities, and 
corrections for branch cuts in the range of integration. 


e Equation solving has been improved. Among these improvements is more versatile 
handling of symbolic matrices in solving systems of linear equations. 


e Interval arithmetic has been improved and extended. 


e Differential equation solving has been enhanced. Symbolic solutions can now be found 
for systems of nonlinear ordinary differential equations, additional Ricatti equations, 
and first- order partial differential equations. 


5.4.2 Package Enhancements 


e There are 11 new packages in such areas as three-dimensional contour plotting, spline 
fitting, variational methods, elliptic integrals, and music and sound. 


e Improvements have been made to existing packages, including graphics, statistics, nu- 
merical limits, and Fourier and Laplace transforms. 


5.4.3 MathLink(R) 


e A new “local” link protocol effectively provides a Mathematica expression type for C 
programmers. 
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New C functions have been added for transferring expressions between links, getting 
low-level device information, storing user data within links, marking positions in in- 
coming data streams, duplicating links, and putting and getting arrays of numbers. 


MathLink now supports the AppleTalk Data Stream Protocol. 


5.4.4 Microsoft Windows Version 


New commands allow the user to find, or find and replace, text; merge and divide cells; 
and call up dialog boxes for editing notebook and default cell styles. 


New commands toggle the notebook window’s ruler and margin markers, the toolbar 
buttons, and the status bar. 


The Action Settings dialog box now provides options for placing print output as it is 
generated, and for placing print output in a separate cell. 


The Graphics Settings dialog box now provides options for specifying the height and 
width of the default graphics size in pixels. 


This version has been extensively tested to enhance efficiency and stability. 


5.4.5 Apple Macintosh Version 


The kernel and the notebook front end are now separate programs communicating via 
MathLink. This allows connecting a front end on one Macintosh to one or more kernels 
on other Macintosh computers; it allows running the kernel without a front end; and it 
results in better handling of kernel interrupts. 


A Function Browser has been added that displays hierarchically organized lists of Math- 
ematica functions, explains these functions, and allows the user to edit function tem- 
plates and paste them into notebook cells. 


Options have been added for printing PostScript renderings of bitmap images, display- 
ing messages, and customizing In/Out prompts and their placement. 


The Graphics Preferences dialog box now provides control of the PostScript image 
operator in the output of DensityGraphics and Raster objects. 


Messages from the front end and kernel can now be displayed in a Messages window 
separate from the notebook window. 


The current time, and the kernel name, can be appended to In/Out prompts. In/Out 
names can be assigned to all “side effects” of evaluations, e.g., information text, graph- 
ics, and messages. An In/Out prompt can be assigned as the keyword of its cell. 


5.4.6 Unix Versions 


MathBook is available under the X Window System; it provides a hypertext interface 
for on-line access to documentation such as The Mathematica Reference Guide and 
Mathematica Warning Messages. 


The X-based PostScript interpreters now support Display PostScript. 
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5.4.7 NeXT Version 


e The kernel and the notebook front end now communicate through MathLink. 


e A Function Browser has been added that displays hierarchically organized lists of Math- 
ematica functions, explains these functions, and allows the user to edit function tem- 
plates and paste them into notebook cells. 


5.5 Maple V Release 2 for PC-compatible computers 


Waterloo Maple Software and the CAN Expertise Centre are proud to announce the availabil- 
ity of a long awaited version of Maple, viz., Maple V, Release 2 on PC-compatible computers 
running MS-DOS or MS-Windows. Release 2 represents a major step forward in all capabil- 
ities of Maple. 

Release 2 incorporates a completely new user interface and many new graphical capabil- 
ities. The new interface under Windows supports worksheets, a help browser, and outputs 
formulae in mathematical notation. High quality hardcopy output for worksheets is available 
in PostScript. 

The new graphical capabilities include contour plots, implicit plots (2D and 3D), vector 
plots (2D and 3D), animation (2D and 3D), and basic support for text, polygons, color, 
lighting, etc. 

As with previous new versions of Maple, Maple V Release 2 also includes many new and 
improved mathematical capabilities. A focus of this version has been given to fixing bugs, 
fixing design problems, and to improve the numerical capabilities, in particular, complex 
numeric arithmetic, numerical linear algebra, and numerical solutions to polynomials. 


SPECAL INTRODUCTION OFFER 


Maple V Release 2 for PC-compatible computers is now available for the special introduction 
price of Dfl. 985,-, which means a 50% DISCOUNT. This offer holds till the end of June 
1993. 


Inquiries about Release 2 can be addressed to 


CAN Expertise Centre 

Kruislaan 413 

1098 SJ Amsterdam 

The Netherlands 

Phone: +31-20-5926050 FAX: +31-20-5924199 


5.6 Derive XM Version 2.5 


Soft Warehouse, Inc. released DERIVE XM, a new Derive product for power users. Derive 
XM requires a 386 or 486 based PC with 2MB of memory. It operates identically to the 
regular Derive and can make use of up to 4 gigabytes of extended memory to solve much 
larger problems. It is distributed on a 3.5” 720K disk with a Derive Manual (available in 
English, German, and Dutch). Derive XM and Derive are sold as different products. 

Derive and Derive XM version 2.55 now got graphics screen capture support built-in! 
Derive and Derive XM now include commands for printing formulas and/or plots to printers 
and for saving a graphics image of the screen to TIF (tagged image format) files. 
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5.7 Nederlandse versie van Derive 


Derive en Derive XM worden sinds 1993 gedistribueerd in een aan de Nederlandse taal 
aangepaste versie: de Derive handleiding, het hulpsysteem en de software-boodschappen 
zijn op initiatief van CAN in het Nederlands vertaald. Hiermee is een mogelijke taalbariére 
overwonnen en is de software gemakkelijker inzetbaar geowrden in het onderwijs. 


5.8 GAP 3.2 


In March 1993, GAP (Groups, Algorithms, and Programming) version 3, Release 2 became 
available. In the previous newsletter, the previous release was described. An updated 
introduction can be obtained from the CAN ftp-archive through the file pub/gap/gap- 
3.2/README. Here, we mention only what’s new in GAP 3.2, which has grown about 
50%. 


e Together with GAP 3.2 are now distributed several *share library packages*. Such 
packages have been contributed by other authors, but the copyright remains with the 
author. Currently there are three packages in the share library. The *ANU PQ* 
package, written by E. O’Brien, consists of a C program implementing a < p >-quotient 
and a < p >-group generation algorithm and functions to interface this program with 
GAP (or Cayley). The *NQ* package, written by W. Nickel, consists of a C program 
implementing an algorithm to compute nilpotent quotients of finitely presented groups 
and a function to call this program from GAP. The *Weyl* package, written by M. 
Geck, contains functions to compute with finite Weyl groups, associated (Iwahori-) 
Hecke algebras, and their representations. 


e It is now possible to extract several elements from a list with a construct similar to 
the one used to extract single elements. This also works recursively, so that it is for 
example possible to extract a submatrix of a matrix. It is also possible to assign several 
elements to a list at once. 


Permutations can now operate on more than 65536 points. 


e Ranges can now also have increments other than 1, ie., a range is now a dense list 
of integers such that the difference between any two consecutive elements is a nonzero 
constant. 


e Strings are now also lists, namely lists of characters, which are a new builtin datatype. 
This makes functions easier to write that deal extensively with strings, such as ’Dis- 
playCharTable’. 


e GAP now supports univariate polynomials over arbitrary coefficient rings. Since the 
coefficient ring may itself be a polynomial ring it is possible to create multivariate 
polynomial rings, though this is not very efficient. Polynomials are implemented in 
the GAP programming language, but there are supporting kernel functions to improve 
efficiency. 


Previously the entries of a matrix had to be among the built-in datatypes, i.e., rationals, 
cyclotomics, and finite field elements. This restriction has been removed, so that it is 
now possible for example to compute with matrices whose entries are polynomials. 


e There is now an implementation of the Dixon-Schneider algorithm, which computes the 
character table of an arbitrary group. 
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e For permutation groups there are new functions to test if a permutation group is solv- 
able, and if so to find a power-commutator presentation. Also there is a new function 
to compute the composition series of a permutation group. 


e The functions to compute presentations for subgroups of finitely presented groups and 
to simplify them are new. 


e There are new functions that work with table of marks, which give a compact description 
of the subgroup lattice of a group. For example there is a function that computes the 
value of the Moebius function for the subgroup lattice of a group with a given table of 
marks. 


e E. O’Brien and C. Rhodes provided a library of 3-groups of size dividing 729. The 
character table library has been extended by about 60 new ordinary tables and about 
200 new modular tables. There is also a data library that contains table of marks for 
various groups, e.g., McL. 


5.9 SIMATH 3.6.1 


In April 1993, the new SIMATH version 3.6.1 became available. SIMATH is a computer 
algebra system for algebraic number theory computations developed by the SIMATH group 
of Prof. Dr. H.G. Zimmer at at the University of Saarland in Saarbriicken, Germany. It is 
written in C, contains an interactive calculator (simcalc) and many C-functions over algebraic 
structures as arbitrary long integers, rational numbers, floating point numbers, polynomials, 
Galoisfields, matrices, elliptic curves, algbraic number fields, modular integers etc. Version 
3.6.1 contains a handbook written in English language. 

SIMATH is available by anonymous ftp from ftp.math.uni-sb.de, ftp.math.orst.edu, and 
canb.can.nl. (you have to get one extra file by an Email-Server, see the installation hints.) 

It runs on Sun 3/140 under SunOS 3.4; Sun 3/80 under SunOS 4.0.3; Sun SPARCstation 
under SunOS 4.1.1; Apollo DN 3000 under SR10.1; Apollo DN 4500 under SR10.1; Apollo 
DN 10000 under SR10.3; Siemens MX2 under SINIX V2.0; Siemens MX2+ under SINIX 
V5.2, Siemens MX300 under SINIX V5.2. and (hopefully) other 32-bit UNIX machines. 


5.10 SymbMath 2.2 


by Dr. Weiguang HUANG 
Dept. Analytical Chemistry, University of New South Wales, 
Kensington, Sydney, NSW 2033, Australia 
Phone: 61-2-697-4643 Fax: 61-2-662-2835 E-mail: w.huang@unsw.edu.au 


SymbMath (an abbreviation for Symbolic Mathematics) is a symbolic calculator that can 
solve symbolical math problems. SymbMath also perform exact numeric computation. It 
can manipulate complicated formulas and return answers in terms of symbols, formulas and 
exact numbers. 

SymbMath is an expert system that is able to learn from users. If users only input one 
formula without writing any code, it will automatically learn many problems related to this 
formula (e.g. it learns the integrals of an unknown function from a derivative of that unknown 
function). 

SymbMath is a programming language in which you can write programs, functions and pro- 
cedures. 
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SymbMath is a computing environment where you can set up, run and document your cal- 
culation. 

SymbMath is a text editor in which you can edit many documents in many windows. 

It runs on small IBM PCs (8086) with 420 KB free memory under MS-DOS. 

Its capabilities include facilities to provide analytical and numerical answers for: 


o Differentiation: regular or higher order, partial or total, mixed and even implicit differ- 
entiation, one-sided derivatives. 


o Integration: indefinite or definite integration, multiple integration, infinity as a bound, 
parametric integration, iterated integration, line or surface integrals, discontinuous or 
implicit integration. 


o Solution of equations: roots of a polynomial, systems of algebraic or differential equations. 


o Manipulation of expressions: simplification, factoring or expansion, substitution, evalua- 
tion, built-in standard or user-defined functions. 


o Calculation: exact and floating-point numerical computation of integer, rational, real and 
complex numbers in the range from minus to plus infinity, even with different units. 


o Limits: real or complex limits, one-sided limits, indeterminate forms. 
o Complex: calculation, functions, derivatives, integration. 

o Sum and product: finite or infinite, partial. 

o Others: series, lists, arrays, vectors, matrices, tables, etc. 

Also included are: 

o Pull-down and pop-up menus, resizable and colourable windows. 
o Twin screen text editor with syntax check. 

o Programming in integrated edit-interpret-debug environment. 

o Reading and generating BASIC and FORTRAN codes. 

o Displaying two-dimensional math notation. 

o On-line help, and on-line manual. 


o Procedural, conditional, iterational, recursive, functional, rule-based, logic, and pattern- 
matching programming. 


o Editable library in ”formula format” or SymbMath language. 
o Learning from user. 
o Plotting functions and data. 


Symbath has three versions: Shareware, Student, and Advanced. The shareware ver- 
sion is available from SIMTEL20 archives on anonymous FTP sites (wsmr-simtel20.army.mil 
as /calculator/sm22a.zip, rana.cc.deakin.oz.au as /huang/sm22a.zip, etc), canb.can.nl as 
/pub/symbmath/sm22a.zip, etc) or by e-mail from listserv@vm1.nodak.edu, or 
listserv@ndsuvm1.bitnet. 


6 | Book Reviews 


6.1 Algorithms for Computer Algebra 


K.O. Geddes, S.R. Czapor & G. Labahn, Algorithms for Computer Algebra, Kluwer Academic 
Publishers, Boston, 1992, 585 p., Dfl. 275, ISBN 0-7923-9259-0 


The Computer Algebra community has been waiting for years for this book to appear. Even 
though one must admit that the wait has been worthwhile, it is not clear that it really should 
have been that long. Nevertheless, the book is a masterpiece and can be recommended to 
everyone interested in the algorithms for computer algebra, either as a reference for further 
research or just to give the casual user an idea why things work as good (or as bad) as they do 
in computer algebra packages. There is one obvious proviso to make in this recommendation: 
the price. The publisher clearly underestimates the market for a book like this, and as usual 
this estimate will be self-fulfilling. 

Let me briefly sketch the contents of the book. 

Chapter 1 (Introduction to computer algebra) has sections on symbolic versus numeric 
computation, history, and as an example a short description of Maple. 

Chapter 2 (Algebra of Polynomials, Rational Functions, and Power Series) gives an in- 
troduction the the mathematical tools. It defines rings, integral domains, UFD’s, Euclidean 
domains and fields. The Euclidean algorithm is treated and its extension. Next are univariate 
and multivariate polynomial domains, the primitive Euclidean algorithm (avoiding divisors). 
Last are rational functions and (extended) power series, and the relation among the different 


domains. 
[ The algorithms are given explicitly 


| in boxes using pseudocode. 

Chapter 3 (Normal Forms and Algebraic Representations) gives an introduction to the 
problems of simplification and uniqueness. E.g. 210° — y190 is simpler in its expanded 
form as given here than in is factored form, even though the latter is more simple from the 
mathematical point of view. Then one turns to the zero equivalence problem, with as an 
example the expression 








sin x 


log(tan(x/2) + see(z/2)) — aresinh( =~ — 


) -l<a2<l. 
The question here is how to decide whether this is zero or not. The chapter closes with the 
internal representation of the data structures. 

Chapter 4 (Arithmetic of Polynomials, Rational Functions, and Power Series) addresses 
the problems encountered in handling large formulae. It gives algorithms for the fast multi- 
plication of polynomials, the Fast Fourier Transform and Newton iteration for power series 
division. 

Chapter 5 (Homomorphisms and Chinese Remainder Algorithms) discusses the possibility 
of avoiding expression swell in integer arithmetic by doing the computations on finite fields 
and lifting the results back to the integer situation. 
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Chapter 6 (Newton’s Iteration and the Hensel Construction) applies the Newton method 
(well known from numerical analysis) to integer problems, a somewhat unexpected field of 
application. Having read this I was able to apply it to my own research, again in a completely 
different direction. This may illustrate that the book by its use of good examples is readable 
for the user of computer algebra who is not an expert in algebra. 

Chapter 7 (Polynomial GCD Computation) applies the techniques developed thus far on 
the GCD problem for polynomials. Needless to say that by now things are getting somewhat 
technical and may not be very interesting for the casual user. 

Chapter 8 (Polynomial Factorization) treats several cases, namely factorization of polyno- 
mials over GF(p™), p prime, and over algebraic extensions (like (Q)(V#) ). This fills for the 
case of GF(p) the gap left in the preceding chapter, where the problem was reduced to factor- 
ization modp. Two algorithms are treated here: Berlekamp’s algorithm and distinct-degree 
factorization. 

Chapter 9 (Solving Systems of Equations) turns to the classical problem of solving both 
linear and polynomial equations. It gives algorithms using resultants to do this. This is much 
refined in 

Chapter 10 (Grébner Bases for Polynomial Ideals) where Grdbner Bases and Buchberger’s 
algorithm are introduced, with applications to Geometry Theorem Proving and Polynomial 
GCS Computation. 

Chapter 11 (Integration of Rational Functions) treats the basic concepts of differential 
algebra and several methods to integrate rational functions. This culminates in 

Chapter 12 (The Risch Algorithm) which gives us an algorithm deciding whether one 
can compute an integral in terms of elementary functions, solving an age-old mathematical 
problem. This has been implemented in several computer algebra packages (AXIOM, Maple). 

If the price was right, the recommendation would be clear: Buy this book! As it stands 
now my only advise is to stay away. from this book, because you might be tempted to buy it 
anyway (at least I was). 

Jan Sanders Amsterdam 


6.2 Ideals, Varieties and Algorithms: An Introduction to 
Computational Algebraic Geometry and Commutative 
Algebra 


David Cox, John Little, and Donal O’Shea, Ideals, Varieties and Algorithms: An Introduc- 
tion to Computational Algebraic Geometry and Commutative Algebra, Undergraduate Text 
in Mathematics, Springer-Verlag, 1992, 513 pages, ISBN 0-387-97847-X and 3-540-97847-X. 


This book is an introduction to algebraic geometry and commutative algebra with emphasis 
on algorithmic aspects. Already after the first introductory chapter, the theory of Grobner 
bases is treated. This forms the basis for almost all algorithms in the rest of the book. Not 
only classical subjects like elimination theory, ideal theory, theory of algebraic varieties, and 
invariant theory of finite groups are dealt with, but also non-classical topics like robotics and 
automatic geometric theorem proving are chosen. In an appendix, the autors discuss the use 
of the computer algebra systems in conjunction with the text. 

It is well-written textbook with many examples that one misses in other books because 
they can only be worked out with a computer algebra system at hand. So it supplements and 
enriches standard introductions of algebraic geometry; certainly, a good choice for students 
and researchers. 

André Heck Amsterdam 
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6.3 Algebraic computing with REDUCE 


M.A.H. MacCALLUM and F.J. WRIGHT, Algebraic computing with REDUCE. Oxford: 
Clarendon Press, 1991. 294 p., £15,-. ISBN 0-19-853444-2, 


Dit boek bevat het eerste deel van het lesmateriaal, dat gebruikt is bij de First Brazilian 
School on Computer Algebra. Het tweede en derde deel dragen resp. als titel Algebraic com- 
puting in general relativety en Tackling mathematical problems with MAPLE. Deze delen zijn 
bij dezelfde uitgever verschenen. 


De eerste ideeén over computeralgebra werden inmiddels 40 jaar geleden gepubliceerd. Sinds- 
dien zijn veel Computer Algebra Systemen ontstaan en ook weer verdwenen. De nu veel 
gebruikte en in ieder geval veel genoemde systemen zijn naast REDUCE b.v. Maple, Mathe- 
matica, Derive, Axiom en Mascyma. Deze opsomming is niet uitputtend en omvat alleen 
CAS die bedoeld zijn voor algemeen gebruik. Naast opmerkelijke onderlinge verschilen 


hebben Axiom, Derive, Macsyma en REDUCE gemeen, dat de codering in een LISP(- 


dialect) gedaan is. Maple en Mathematica zijn in C gecodeerd. De op LISP gebaseerde 
systemen zijn (veel) eerder ontstaan dan het C tweetal, waarvan distributie op duidelijk com- 
merciéle leest geschoeid is. Voor REDUCE geldt dat nog steeds niet. Het is vrijwel gelijk 
met LISP ontstaan, aanvankelijk als eenmans-product. De eerste versie van REDUCE werd 
door A.C. Hearn, die nog steeds de eindverantwoording voor het systeem draagt, ontworpen 
voor berekeningen in de hoge energie fysica. Het was toen een concurrent van Veltman’s 
SCHOONSCHIP, nu geleidelijk aan opgevolgd door Vermaseren’s FORM. Hearn werkt nu 
voor Rand, een non-profit organisatie, en distribueert vandaar uit REDUCE, samen met an- 
deren, voor een kostendekkend bedrag. De ontwikkeling, die het systeem in de loop der jaren 
doorgemaakt heeft zou beschreven kunnen worden als een codperatieve onderneming van ve- 
len. Van de faciliteiten van de huidige versie van REDUCE is het grootste deel door Europese 
onderzoekers gemaakt en ter beschikking gesteld. Dit was vooral mogelijk omdat REDUCE 
een open systeem is; de tape bevat naast de source code instructies voor installatie. Dit 
open karakter maakt het mogelijk het systeem naar believen te veranderen en uit te breiden. 
Hiervoor is natuurlijk wel kennis nodig over de structuur en samenhang tussen de samenstel- 
lende delen. Een andere oorzaak -en tevens pluspunt- is de reeds lang bestaande portabiliteit 
van REDUCE. Hearn en zijn naaste medewerkers hadden daartoe een eigen LISP-standaard 
ontwikkeld. Dit wat magere LISP werd als basis gebruikt voor hun PSL, Portable Standard 
LISP. Dit PSL tenslotte is gebruikt om de communicatievoorzieningen tussen gebruiker en 
systeem vorm te geven. REDUCE kent twee gebruikslagen. Het hoogste niveau, de Alge- 
braic Mode, laat interactie met de machine toe via commando’s, waarvan de syntactische 
constructies sterk lijken op wat gangbaar is bij programmering met Pascal. Een niveau lager, 
de Symbolic Mode, laat LISP-achtige programmering toe. De opdrachten kennen ook hier 
een Pascal-verpakking. Bij dit niveau wordt veelal gesproken van RLISP. RLISP wordt i.h.a. 
gebruikt voor systeem-ontwikkeling. 


Een CAS kan opgevat worden als een intelligente calculator. Dan is de omvang van en toe- 
gang tot de in het systeem opgeslagen wiskundige kennis het toetsingscriterium. Een CAS 
kan ook gezien worden als een ontwikkelomgeving. Het maken en implementeren van nieuwe 
mogelijkheden staat dan centraal. Bij de beoordeling van een boek, dat een CAS als onder- 
werp heeft moeten beide aspecten aandacht krijgen. Bij een open systeem als REDUCE is 
met name informatie over de zorgvuldigheid die een gebruiker moet hanteren om het systeem 
niet in wanorde te brengen van groot belang. Verder dient bezien te worden of de verstrekte 
informatie iets toevoegt aan bestaand materiaal, hetzij uit didaktische, hetzij uit louter in- 
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formatieve overwegingen. De REDUCE tape bevat een handleiding. Deze vooronderstelt 
feitelijk dat de gebruiker LISP al kent en zelf in staat is te begrijpen hoe RLISP gemaakt en 
toegevoegd is. Dat betekent dat nieuwe gebruikers steeds eerst die barriére moeten nemen. 
De conclusie moge duidelijk zijn: een boek over REDUCE is pas echt interessant als we ook 
didaktisch verantwoorde informatie over RLISP en het gebruik ervan aantreffen. 


Algebraic computing with REDUCE is opgebouwd uit negen hoofdstukken en een appendix 
met enige technische details. Ieder hoofdstuk wordt afgesloten met een aantrekkelijke verza- 
meling oefeningen, die van makkelijk oplopen tot vrij pittig. De tekst is gelardeerd met vele 
voorbeelden en praktische tips. 

In hoofstuk 1 wordt REDUCE gepresenteerd als calculator. De lezer wordt aan de hand 
meegenomen, ervan uitgaande dat hij over een bedrijfsklare versie beschikt. Stap voor stap 
wordt kennisgemaakt met het vraag-antwoord mechanisme, waar nodig voorzien van achter- 
grond informatie, b.v. over getalrepresentatie en arithmetiek, differentiatie en integratie, 
polynomen en rationale funkties. 

Hoofdstuk 2 handelt over een wat genuanceerdere en kritischer commando-verwerking. 
Een REDUCE-opdracht ondergaat een standaardbehandeling: de invoer wordt ”vereen- 
voudigd”, opgeslagen en het resultaat wordt ih.a. weer zichtbaar gemaakt. Deze vereen- 
voudiging van aangeboden wiskundige uitdrukkingen wordt op de invoer steeds volgens 
dezelfde regels uitgevoerd. Deze regels kunnen door de gebruiker aangevuld en bijgesteld wor- 
den. Deze, toch wat paternalistische aanpak -wij weten wat goed is- is in nieuwere systemen, 
als b.v. Maple, ten dele losgelaten. Hier kan bij iedere nieuwe opdracht, via funktie-applicatie 
en parameterkeuze, een locale verwerkingsomgeving gemaakt worden. Een trend, die in de 
nieuwste versie van REDUCE ook in gang gezet is. 

Hoofstuk 3 wordt gebruikt om een overzicht te geven van de diverse taalconstructies, die 
REDUCE, opgevat als programmeertaal, kent. In hoofdstuk 4 wordt hierop voortgebouwd bij 
de behandeling van procedures. Procedures kunnen in een interactief systeem als REDUCE 
gebruikt worden om de omgeving te verrijken met samengestelde, soms complexe, behandel- 
ingsvoorschriften. Mits syntactisch correct kunnen deze na invoering toegepast worden bij 
vormgeving van volgende opdrachten. Ook operatoren, initiéel niet-gespecificeerde funkties, 
kunnen hierbij een rol spelen en worden daarom ook in dit hoofdstuk behandeld. 

Hoofdstuk 5 wordt gebruikt om de blik naar binnen te richten: de relatie tusen (R)LISP 
en REDUCE komt aan de orde. Hierop wordt in de hoofdstukken 6 en 9 voortgebouwd, bij 
de behandeling van RLISP. 

Bij de hoofdstukken 7 en 8 tenslotte komt REDUCE weer als calculator aan bod. Hoofd- 

stuk 7 staat wat los van de rest van het boek. Polynoomfactorisatie en integratie worden hier 
van enige context voorzien. Deze onderwerpen vormen aansprekende resultaten van door het 
ontstaan van CAS geinspireerd onderzoek. Hoofdstuk 8 is gewijd aan door gebruikers van 
REDUCE ontwikkelde en ter beschikking gestelde applicatie-pakketten, zoals GENTRAN, 
t.b.v. code-generatie, en het Grdbner-pakket. Bij de in dit boek behandelde 3.3 versie van 
REDUCE konden deze pakketten via electronische postvoorzieningen verkregen worden. Bij 
de nieuwe 3.4 versie is een 15-tal van deze pakketten al meegeleverd op de systeem-tape. 
De postvoorziening wordt nu gebruikt voor communicatie over en t.b.v. gebruik van meer 
recente of nog experimentele voorzieningen. 
Het boek is goed geschreven en kan uitstekend voor zelfstudie gebruikt worden. De opgaven 
kunnen goed dienst doen als aanvulling op lesmateriaal. Het is bij mijn weten het enige 
engelstalige boek dat nuttige en uitvoerige informatie over RLISP bevat. Er is eigenlijk 
maar €én ding jammer: het boek verscheen kort voordat de nieuwe versie van REDUCE ter 
beschikking kwam. Maar bij herdruk zal dit euvel zeker verholpen worden. Kortom: van 
harte aangebovelen. J.A. van Hulzen, Enschede. 


7 | Book Announcements 


7.1 Computeralgebra met Maple 


Metha Kamminga van Hulsen, Computeralgebra met Maple, Academic Service, 1993, 187 p., 
Dfi. 38,—, ISBN 90 6233 923 9. 


Gezien de ontwikkelingen van de laatste jaren op het gebied van computeralgebrasystemen, 
ligt het voor de hand dat elke student, HBO of universitair, die te maken heeft met wiskunde, 
tenminste kennis moet maken met een dergelijk systeem in zijn of haar opleiding. 

Dit boek is tot stand gekomen bij de afdeling werktuigbouwkunde van de Noordelijke 
Hogeschool Leeuwarden, waar in 1991 computeralgebra in de wiskundeles is ingevoerd. 

Er is gekozen voor MAPLE, een volwaardig computeralgebrasysteem, waarin ook gepro- 
grammeerd kan worden. Dit systeem is ook al bij bedrijven bekend en heeft vele toepassingen 
om analytisch te rekenen. Denk aan formele algebra, formulemanipulatie, differentiren en 
integreren, oplossen van vergelijkingen en stelsels vergelijkingen en het oplossen van differen- 
tiaalvergelijkingen. 

Deze handleiding kan parallel aan de wiskundelessen gebruikt worden op HBO en univer- 
sitair niveau. De hoofdstukken kunnen zelfstandig doorgewerkt worden en de opgaven aan het 
eind van ieder hoofdstuk zijn bedoeld om met het computeralgebrasysteem te oefenen. (Er 
wordt bij dezelfde uitgever een aparte bundel uitgegeven met uitgebreide vraagstukken die bij 
voorkeur met een computeralgebrasysteem aangepakt moeten worden). Na het doorwerken 
van deze handleiding is men in staat om goed met het handboek van MAPLE om te gaan en 
zich verder in het gebruik van computeralgebrasystemen in het algemeen te bekwamen. 

Veel onderwerpen die MAPLE ook aankan komen in dit boek niet aan de orde zoals de 
pakketten: statistiek, lineaire optimalisering, groepentheorie, Euclidische meetkunde, projec- 
tieve meetkunde enzovoort. Bovendien valt het programmeren binnen MAPLE ook buiten 
het bestek van dit boek. 

Voor gevorderden in de wiskunde is deze handleiding ook snel door te werken en is zelfs 
leesbaar zonder een computer bij de hand. De in en uitvoer is zo afgedrukt zoals deze ook te 
zien is op het scherm. 

Dit boek behandelt beslist geen theorie van de wiskunde maar het sluit wel aan bij elk 
willekeurig theorieboek op HBO of universitair niveau. 

Korte beschrijving van de inhoud: 


e Hoofdstuk 1 en Hoofdstuk 2 gaan over basisvaardigheden. Zoals de helpfunctie, enkele 
commando’s voor formulemanipulatie, vergelijkingen oplossen, limieten, sommatie en 
differentiéren. 


e Hoofdstuk 3 en Hoofdstuk 4 gaan over plotten en integreren: Onbepaald integreren, de 
bepaalde integraal, oneigenlijke integralen en numerieke benaderingen. 


e Hoofdstuk 5 is een uitbreiding van het differentifen. Hierin komen dingen aan de 
orde zoals de hogere afgeleiden en partiéle afgeleiden, reeksontwikkeling, relaties in 
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parametervorm, impliciete functies. 


Hoofdstuk 6 is een hoofdstuk waarin een aantal tips en trucs worden besproken die bij 
het doorwerken van de eeste vijf hoofdstukken ook al van nut kunnen zijn. 


In Hoofdstuk 7 worden de complexe getallen behandeld. 


Hoofdstuk 8 is een overzicht van wat er allemaal mogelijk is bij de formulemanuipulatie. 


In Hoofdstuk 9 en 10 wordt lineaire algebra behandeld. Hier komen de bewerkingen 
met vectoren en matrices aan de orde. Determinanten, eigenwaarden en eigenvectoren, 
lengte, hoek en in en uitproduct. 


e Hoofdstuk 11 gaat over differentiaalvergelijkingen. 


In Hoofdstuk 12 wordt de Laplacetransformatie behandeld met alles wat daarbij hoort, 
zoals de Diracpuls, de Sprongfunctie en het oplossen van differentiaalvergelijkingen met 
Laplace. 


Tot slot komen in Hoofdstuk 13 de meervoudige integralen aan de orde. Ook wordt op 
de poolcordinaten ingegaan met de daarbijbehorende grafieken en nog een uitbreiding 
van de plotfunctie met cilindercordinaten en bolcordinaten. 


7.2 Introduction to Maple 


André Heck, Introduction to Maple, Springer Verlag, 1993, 497 p., ISBN 0-387-97662-0 and 
3-540-97662-0. 


‘This long awaited book will be available in June 1993. It is the first real introduction to 
Maple V, Release 2 and provides a gentle presentation of this computer algebra system. It 
intends to teach not only what can be done by the system, but also how and why it can 
be done. Emphasis is more on understanding Maple than on factual knowledge of built-in 
facilities. To this end, the book contains many examples and exercises, both elementary and 
more sophisticated. They stimulate the reader to use Maple and encourage to find one’s way 
through the system. 


In this book, the usage of Maple as a programming language is not discussed at a higher 
level than that of defining simple procedures and using simple language constructs. However, 
the Maple data structures are discussed in great detail because good understanding of them 
is necessary for manipulating and simplifying expressions effectively. This also forms a good 
starting point to acquaint oneself further with Maple as a programming language. 


Contents 


1. Introduction to Computer Algebra 

What is Computer Algebra? —- Computer Algebra Systems — Some Properties of Computer 
Algebra Systems — Advantages of Computer Algebra — Limitations of Computer Algebra 
— Maple. 

2. The First Steps: Calculus on Numbers 

Getting Started — Getting Help — Integers and Rational Numbers — Irrational Numbers 
and Floating-Point Numbers — Algebraic Numbers — Complex Numbers — Exercises. 

3. Variables and Names 

Assignment and Evaluation — Unassignment — Full Evaluation — Names of Variables — 
Basic Data Types — Exercises. 
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4. Getting Around with Maple 
Input and Output — The Maple Library — Reading and Writing Files — Formatted I/O — 
Code Generation — Changing Maple to your own Taste — Exercises. 


5. Polynomials and Rational Functions 
Univariate Polynomials — Multivariate Polynomials — Rational Functions — Conversions 
— Exercises. 


6. Internal Data Representation and Substitution 
Internal Representation of Polynomials — Generalized Rational Expressions — Substitution 
— Exercises. 


7. Manipulation of Polynomials and Rational Expressions 
Expansion — Factorization — Canonical Form and Normal Form — Normalization — Col- 
lection — Sorting — Exercises. 


8. Functions 
Mathematical Functions — Arrow Operators — Maple Procedures — Recursive Procedure 
Definitions — unapply — Operations on Functions — Anonymous Functions — Exercises. 


9. Differentiation 
Symbolic Differentiation — Automatic Differentiation — Exercises. 


10. Integration and Summation 
Indefinite Integration — Definite Integration — Numerical Integration — Integral Transforms 
— Assisting Maple’s Integrator — Summation — Exercises. 


11. Truncated Series Expansions, Power Series, and Limits 
Truncated Series Expansions — Power Series — Limits — Exercises. 


12. Composite Data Types 
Sequence — Set — List — Array — convert and map — Exercises. 


13. Simplification 
Automatic Simplification — expand — combine — simplify — convert — Trigonometric 
Simplification— Simplification w.r.t. Side Relations — Exercises. 


14. Graphics 

Some Basic Two-Dimensional Plots — Options of plot — The Structure of Two-Dimensional 
Graphics — Special Two-Dimensional Plots — Plot Aliasing — A Common Mistake — Op- 
tions of plot3d — The Structure of Three-Dimensional Graphics — Special Three-Dimensional 
Plots — Animation — Exercises. 


15. Solving Equations 

Equations in One Unknown — Abbreviations in solve — Some Difficulties — Systems of 
Equations — The Grobner Basis Method — Numerical Solvers — Other Solvers in Maple — 
Exercises. 


16. Differential Equations 
First Glance at ODEs — Analytic Solutions — Taylor Series Method — Power Series Method 
—~ Numerical Solutions — Perturbation Methods — Liesymm — Exercises. 


17. Linear Algebra: Basics 
Basic Operations on Matrices — Last Name Evaluation — The Linear Algebra Package — 
Exercises. , 


18. Linear Algebra: Applications 

Kinematics of the Stanford Manipulator — A 3-Compartment Model of Cadmium Transfer 
~~ Molecular-orbital Hiickel Theory — Prolate Spheroidal Coordinates — Moore-Penrose 
Inverse — Exercises. 
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7.3 Economic and Financial Modeling with Mathematica 


Economic and Financial Modeling with Mathematica, Hal Varian, editor, TELOS/Springer- 
Verlag, 1993, ISBN 0-387-97882-8. (Includes disk.) 


This book consists of 17 chapters describing how Mathematica can be used in economic and 
financial analysis. A DOS 3.5 inch disk contains the Mathematica Packages described in the 
book. 


Contents 

Hal Varian (U Michigan), Symbolic Optimization 

Todd Kaplan & Arijit Mukherji (UMinn), Designing an Incentive Compatible Contract 
John Eckalbar (Cal State, Chico), Economic Dynamics 

Ken Judd & Sy-Ming Guu (Stanford), Perturbation Solution Methods for Economic Growth 
Models 

Asahi Noguchi (Senshu U, Japan) General equilibrium models 

Gary S. Anderson (Maryland), Symbolic Algebra Programming for Analyzing the Long Run 
Dynamics of Economic Models 

John Dickhaut & Todd Kaplan (UMinn), A Program for Finding Nash Equilibria 

Michael Carter (Christchurch), Cooperative Games 

Michael Steele & Bob Stine (UPenn), Mathematica and Diffusions 

Wilfrid Kendall (Warwick), Itovsn3: Doing Stochastic Calculus with Mathematica 

Colin Rose (Sydney) Bounded and Unbounded Stochastic Processes 

Ross Miller (General Electric), Option Valuation 

Steve Brown (NYU) Nonlinear Systems Estimation: Asset Pricing Model Application 
David Belsley (Boston College) Econometrics.m: A Package for Doing Econometrics in 
Mathematica 

Eduardo Ley & Mark Steel (Don Carlos III, Madrid), Bayesian Econometrics: Conjugate 
Analysis and Rejection Sampling 


Bob Stine (UPenn), Time series Models and Mathematica 


Robert Korsan (Decisions! Decisions!), Decision Analytica: An Example of Bayesian Infer- 
ence and Decision Theory Using Mathematica 


7.4 Computer Algebra in Industry 


Computer Algebra in Industry, Proceedings of the 1991 SCAFI Seminar at CWI (Amster- 
dam), Arjeh M. Cohen (ed.), Wiley & Sons, 1993, 262 p., £24.95, ISBN 0-471-93829-7. 


CAN is pleased to announce the publication of the book that is the outgrowth of the pro- 
ceedings of the SCAFI’91 seminar. SCAFI stands for Studies in Computer Algebra For 
Industry, a project of Computer Algebra Nederland (CAN) to show how computer algebra in 
industry can be useful and cost effective. The book provides nine case studies on industrial 
applications of computer algebra, developed in collaboration with a number of European 
companies. They cover important topics which include robotics, fluid mechanics, elasticity 
theory; geodesy, channel transmissions and electronic circuit development. Also included is 
a general introduction to computer algebra and two contributions explaining how numerical 
mathematics software can successfully interact with computer algebra. 


Contents + Abstracts 
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General Introduction 
A.M. Cohen, J.H. Davenport, and A.J.P. Heck, An overview of computer algebra, 1 — 52 


The goal of this paper is to briefly describe what computer algebra is about, present a little 
history, some examples of computer algebra use, and discuss some new directions of this new 
and promising field. 

Liquids 

P. Boomkamp and R. Miesen, Systems for doing mathematics by computer: convenient tools 
in luid mechanics, 53 — 70 


Systems for doing mathematics by computer are convenient tools for mathematical modelling 
in general, and in doing fluid mechanical analyses in particular; they enable us to handle 
problems that would be too complicated or too expensive without them. In the present 
paper, this thesis is illustrated by the study of a fluid mechanical problem that is related to 
the transport of very viscous oil. The essence of the problem is a stability analysis of the 
interface of a two-fluid Poisseuille flow in a pipe, i.e., core-annular flow, with a very viscous 
fluid (oil) as the core and a less viscous fluid (water) as the annulus. It is shown that both 
axisymmetric as well as non-axisymmetric interfacial wave modes will be present. 


R.M.J. Schielen, An amplitude equation for the non-linear analysis of river beds, 71 — 86 


An amplitude equation for weak non-linear interactions of waves is derived in order to describe 
the time and spatial evolution of a so-called alternate bar pattern that arises, in auspicious 
conditions, on a river bed. The coefficients of the equation are determined in an asymptotic 
expansion of the drag coefficient, which is small in natural circumstances. This expansion 
has the advantage that the influence of the model parameters on the amplitude equation can 
be seen immediately. 


A.G. Schlijper and S.M. Verduyn Lunel, Higher order interactions in statistical mechanics, 
87 ~ 100 


Wetting and adsorption are industrially relevant phenomena. A project at Koninklijke/Shell 
Exploratie en Produktie Laboratorium (KSEPL) aims at the development of a statistical 
mechanical theory which describes these phenomena at a molecular level. Computer algebra 
will be an invaluable tool to provide higher order approximations to the exact statistical 
mechanical theory. 


Mechanics 
P.J. Cool, Equations of motion of a stabilised platform, 101 - 116 


This paper is about the derivation of the equations of motion for a stabilised radar antenna 
on a platform installed on a ship. Two models for the system will be derived, one for a test 
situation on land and one for the situation on sea. In fact the derivation procedure consists of 
writing down Newton’s Law for each part of the system and of connecting all parts involved. 
Because this procedure is a lot of work, it is done with the assistance of the computer algebra 
package Maple. 


M.J. Gonzalez-Lopez and T. Recio, The ROMIN inverse geometric model and the dynamic 
evaluation method, 117 — 142 


The inverse geometric model of a robot manipulator involves finding the joint angles for a 
given place (position and orientation) of the manipulator tip. This problem can be stated 
from an algebraic point of view in such a way that its solution consists of solving a system of 
polynomial equations with parametric coefficients in which the variables represent the joint 
angles and the parameters the place of the tip. The method of dynamic evaluation applied 
to the solution of polynomial systems with parametric coefficients provides a useful tool to 
attack this problem in practice because we get a finite collection of completely triangular 
systems (also with inequalities) whose union is equivalent to the original one. Several other 
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symbolic computation methods for solving the same problem are discussed. The dynamic 
evaluation method has been applied to solve the inverse geometric model of the manipulator 
ROMIN, developed as part of a joint project between the University of Cantabria, the enter- 
prise ENSA (Santander) and the Instituto de Automatica Industrial of Madrid, to perform 
manipulation tasks in high radioactivity zones of a nuclear plant. We also include a detailed 
description of the algorithm for this particular case, together with some numerical examples 
of its performance. 


P. Rideau, Computer algebra and mechanics. The James software, 143 — 158 


AEROSPATIALE has developed a software package called JAMES which is used for the com- 
putation of dynamic equations of multi-body systems. This software is based on the Maple 
manipulation software. This paper shows an application of this software to the deployment 
of a solar array mockup. 


Electronics 
A.M. Cohen, L.J. van Gastel, and A. Jacobs, Error bursts in transmission channels, 159 — 
174 


A Markov model due to Gilbert and Elliott is used for describing error bursts in transmission 
channels. We compute several quality measures of the CCITT for digital transmission in 
terms of the parameters of the Markov model. Computer algebra turns out to be useful for 
manipulation of large formulae, for Taylor expansion and for handling special functions. 
T.H. Koornwinder, On the recurrence relations connected with Elliott’s formula: Appendiz to 
"Error bursts in transmission channels”, 175 — 178 


A closed form is given for the expressions G(m,n) and B(m,n) of the previous paper. The 
resulting formula however is not suitable for efficient computations, except in some special 
cases, which are discussed at the end of this appendix. 


R.M. Aarts and §.M. Verduyn Lunel, Symbolic analysis of loudspeaker cross over filters, 179 
— 186 


Computer algebra can be used to compute the symbolic transfer function of a loudspeaker 
cross over filter as a function of frequency, filter components and loudspeaker impedance. We 
describe a key ingredient of this kind in a program to analyse RLC-loudspeaker cross-over 
filters developed and implemented at the Philips Research Lab., Eindhoven. 


Miscellaneous 

G. Braat, L. van Gastel, A. Heck, L. Saes, and J. Sanders, Strain in rubber coated rollers, 
187 — 202 

For two steel rollers with several rubber layers that are pressed together, a system of equations 
for the displacement and strain is obtained along the lines of the classical theory of linear 
viscoelasticity. The use of computer algebra for linear algebra and code generation enables 
us to solve the equations, even if the number of layers becomes larger. 

A. Heck, Transformation between geocentric and geodetic coordinates, 203 — 220 


Currently, iterative algorithms are much used in computational geodesy for solving the con- 
version from geocentric to geodetic coordinates. We show that use of computer algebra readily 
leads to a closed form solution, and compare our approach with other solution methods. 


Numeric Interfaces 
M. Dewar, Integrating symbolic and numeric computation, 221 — 232 
There is an increasing awareness of the need to integrate both symbolic and numeric methods 


in problem-solving systems. This paper gives some motivation for this and describes both 
the available technology and what may become available in the future. 


M.G. Richardson, The IRENA User Interface to the NAG Fortran Library, 233 — 244 
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The NAG Library provides powerful numerical analysis routines. IRENA significantly sim- 
plifies the use of these tools. 


7.5 Teaching Mathematics with Derive 


Teaching Mathematics with Derive, Proceedings of the KREMS ’92 conference, J. Bohm 
(ed.), Chartwell-Bratt, 1993, 300 pages, £13.95, ISBN 0-86238-319-6. 


This Proceeding of the KREMS ’92 conference is an illustrated guide packed full of sound 
advice and teaching materials. Top European educationalists suggest ways tuse Derive to 
teach mathematics (particularly to 12-18 years ols). Full of helpful suggestions for education- 
alists, it is thougt-provoking impact on the curriculum. If you are teaching mathematics, at 
school, college, or university, you need to know about this new software. Its full impact on 
mathematics in general and teaching in particular will be tremendous. New approaches are 
explored in this book. 


7.6 Partial Differential Equations with Mathematica 


Dimitri Vvedensky, Partial Differential Equations with Mathematica, Addison-Wesley, 1993, 
465 pages, ISBN 0-201-54409-1 


This book comprises a course on partial differential equations for physicists, engineers, 
and mathematicians. Uses a geometric approach in providing an overview of mathemati- 
cal physics. Uses Mathematica to perform complex algebraic manipulations, display simple 
animations and 3D solutions, and to write programs to solve differential equations. 


7.7 AXIOM, the Scientific Computation System 


Richard D. Jenks, Robert S. Sutor, AXIOM, the Scientific Computation System, Springer 
Verlag, 1992, 742 pages, ISBN 0-387-97855-0. 


This book is the user guide to the AXIOM system, describing it both as an advanced 
symbolic calculator and as a programming language for mathematical problem solving on 
computers. 


7.8 Een computerprogramma snelle Fourier-transformaties 


René Beerends, André Heck, Evert van de Vrie, Een computerprogramma snelle Fourier- 
transformaties, CONCAVO-publikatie 4, Amsterdam, 1993. 


In deze publikate wordt een uitgebreide beschrijving gegeven van de werking van het com- 
puterprogramma Discrete Fourier-transformaties passend bij de Open Universiteitscursus 
Fourier- en Laplace-transformaties. Toepassingsmogelijkheden van het programma worden 
behandeld. Tevens wordt de software vergeleken met de computeralgebra Derive, Maple en 
Mathematica. 
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7.9 Computer Algebra in Deutschland 


Computer Algebra in Deutschland: Bestandsaufnahme, Méglichkeiten, Perspektiven. Fach- 
gruppe Computeralgebra der GI, DMV, GAMM, Passau, 1993 


Auf 320 Seiten haben zu diesem Report 99 Autoren mit Grundsatzlichem, mit Ubersichtsartikeln, 
Anwendungsbeschreibungen, Berichten zu den Systemen und sonstigen Informationen beige- 
tragen. Wir zitieren aus dem Vorwort von Professor Engeler aus Ziirich: 


In dem vorliegenden Report wird meines Wissens zum ersten Mal versucht, einen 
méglichst umfassenden Uberblick iiber dieses Gebiet zu geben. Hierzu gehdren 
neben einer Ubersicht iiber die Forschungsschwerpunkte einige Berichte iiber An- 
wendungen, eine Beschreibung aktueller Computeralgebrasysteme sowie Informa- 
tionen tiber Fachtagungen und Publikationen. Der letzte Teil enthalt schlieSlich 
eine Aufstellung der Aktivitaéten der Fachgruppe Computeralgebra sowie eine 
Liste von CA-Arbeitsgruppen in Deutschland. 


Der Report kann zum Preis von DM 25,00 bei der DLGI, einer Tochter der GI, durch Brief, 
Telefon, Telefax, elektronische Post oder aber durch Uberweisung mit entsprechendem Ver- 
merk auf dem Uberweisungstrager bestellen werden: 


DLGI Dienstleistungsgesellschaft fiir Informatik mbH, 

Godesberger Allee 99, W-5300 Bonn 2 (ab 1.7. 53175 Bonn) 

Telefon 0228-95994-11 oder -16 

Telefax 0228-95994-20 

Elektr. Adr.: gibonn@gmd.de (mit der Bitte um Weiterleitung zur DLGI) 
Kontonummer 30403, Sparkasse Bonn, Bankleitzah] 380 500 00 

Vermerk: Computeralgebra-Report 


7.10 Other CA-Books 


e D.C. Arney, Exploring Calculus with DERIVE, Addison-Wesley, 1992, 166 pages, ISBN 
0-201-52839-8. 


e M. Kofler, Mathematica. Einftihrung und Leitfaden fiir den Praktiker, Addison-Wesley, 
Bonn, 1992. 


e Margret Hoeft, Laboratories for Calculus I using MATHEMATICA, 1992, Addison- 
Wesley, 1992, ISBN 0-201-54345-1. 


e J.K. Finch, M. Lehman, Exploring Calculus with MATHEMATICA, Addison-Wesley, 
1992, ISBN 0-201-55572-7. 


F.W. Hehl, V.Winkelmann and H.Meyer, Computeralgebra: Ein Kompaktkurs tiber die 
Anwendung von REDUCE, Springer-Verlag, 1992. 
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8 | New Features in Maple V, Release 2, 
PART IT 


Symbolic Computation Group} 





8.1 Preface 


This is a continuation of the list of new features of Maple V, Release 2, published in the 
previous newsletter 9, §12. Previously were described: 


e Development of Improved User Interfaces, 
e Graphics, and 


e New and Enhanced Mathematical Facilities. 


8.2 New and Enhanced Integration Facilities 


8.2.1 Dirac and Heaviside 
Integration now knows about the Heaviside and Dirac functions, e.g. 
> int(f£(x)*Dirac(x) ,x=-infinity..infinity) ; 
£(0) 
> int (f(x) *Heaviside(x) ,x=-3..4); 


4 
/ 


| 
| £(x) dx 
| 


/ 
0 
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8.2.2 Elliptic integrals 


Integration now recognizes Elliptic integrals and reduces them to a normal form in terms of 
Legendre functions 
> £ := sqrt( (1-x°2) * (1/25 + x°2/4) ); 


2 1/2 2 1/2 
f£ := 1/10 (1 - x ) (4+ 25x) 


> int(1/f,x=1/2..1); 
10- 1/2 1/2 1/2 
---- 29 LegendreF(1/2 3. , 5/29 29) 
29 
> evalf("); 


2.327355206 


8.2.3 Algebraic functions 


Integration of algebraic functions now uses the Risch-Trager algorithm. This includes a 
full implementation of the algebraic part including integrands which involve parameters and 
algebraic numbers, for example 


> alias (alpha=RootOf(y*3-x*2-a,y); # i.e. alpha = (x*2+a)*(1/3) 
> £ := (4*alpha*2#x°3+ (5*x°4+3*x"2*a) *alpha-3*x*2-3%a) /(x*2+a) /x°2; 


2 3 4 2 2 
4 alpha x + (5x +3x a) alpha-3x -3a 
2 2 
(x +a) x 
> int(f,x); 
x - 1 2 
SG .cse== + 3 x alpha + 3 alpha 
x 


and an implementation of the transcendental part for algebraic functions which do not 
contain parameters, e.g. 


> alias (alpha=RootOf(y*2-x"3-1,y)); # i.e. beta = sqrt(x*3+1) 
> £ := (3#x°4-6*x°3+5%x"2-24x+2) / (x7 5-2#x74-x+2*x°3+1) *beta 


> +(34x°44+5%x°2-5%x~3-2#x) / (x7 5-2*x74-x+24x73+1); 
4 3 2 4 2 3 
(3x -6x +5x -2x+2) beta 3x +5x -5x -2x 
fo tm ---~-----~~ ee + ee we ne 
5 4 3 5 4 3 
x -2x -x+2x +1 X= 22 =x +22 #4 
> int(f,x); 
2 3 
x +x +1-2x beta 
2 beta - 1n(---------------------- ) 
‘ 3 2 2 
(x - x + 1) 


8.2.4 Integration of the W function 


Indefinite integration now uses an inverse function transformation to transform integrals 
involving the W function (where W(2) is the function defined to satisfy eW()W(zx) = x ) of 
linear arguments to elementary functions. Hence the following integrals are now computed 


See eee URE ee 


Ri 
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> int (x/W(x),x); 


2 2 
: x x 
or ee, 
W(x) 2 
W(x) 
> int(1/W(x) ,x); 
x 
---- - Ei(1, - W(x)) 
W(x) 


8.2.5 The exponential integral, error function, and related integrals 


Indefinite integration of the forms g(x) sin(f(z)) and g(x) cos(f(x)) where g(x) is a rational 
function in z and f(x) is a quadratic polynomial in x are expressed in terms of the Fresnel 
integrals FresnelC(x) and FresnelS(x), and the cosine and sine integrals Ci(x) and Si(x). Note 
that, Maple V can do these examples, but the answer involves complex error functions and 
exponential integrals of complex arguments 


> int( 5*x*4*sin(2*x72+1), x ); 


3 2 15 2 
- 5/4 x cos(2x + 1) +’---- x sin(2x + 1) 
16 
15 1/2 / x x \ 
on |cos(1) FresnelS(2 ----- ) + sin(1) FresnelC(2 ----- )| 
32 1/2 1/2 | 
\ Pi Pi i 


> int( sin(u*x*2-v)/x, x ); 
2 2 
1/2 cos(v) Si(u x ) - 1/2 sin(v) Ci(u x ) 


8.2.6 New classes of definite integrals 


Several new classes of infinite indefinite integrals are now handled. These are computed by 
evaluating derivatives of special functions and exploiting symmetry. For example, for p > 0 
and q > 0, the class 


; a—1,—px*—qar-* 
ee dz 
0 


in general can be expressed in terms of the Bessel K function, for example 


> assume(p>0) ; 
> assume (q>0) ; 
> int( x°2*exp(-p*x73-q/x°3), x=0..infinity ); 


1/2 1/2 1/2 
q” BesselK(1, 2 p~ q 3 


Note, the appearance of p~ and q” in this answer indicates that assumptions about p and 
q have been made. Also, several classes involving Bessel functions e.g. for u > 0,v > 0, the 
class 


/ ee 7, (a has 
0 


in general results in the confluent hypergeomic function, for example 
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> int( exp(-t)*t°2*BesselJ(1,sqrt(t)), t=O..infinity ); 


15 1/2 
---- Pi hypergeom([7/2], [2], -1/4) 
16 


Also classes of functions involving the error function can now be handled. For example, for 
u>O0 


> assume(u>0); 
> Int( exp(-u*x*2)*erf(b*x), x=O0..infinity ) = 
> int( exp(-u*x*2)*erf(b*x), x=0..infinity ); 


b 
infinity arctan (----~ ) 
1/2 
| 2 u 
| exp(= n° x } arith x) de = <=<-++----4=- 
| 1/2 1/2 
‘d Pi um 
0 


8.2.7 Handling singularities in an indefinite integral 


Definite integration now uses a new algorithm to resolve integrals which have discontinuities 
(removable or non-removable) in the result from indefinite integration. For example, 
> int( 1/(5+3*cos(x)), x=0..13*Pi/2 ); 


1/2 arctan(1/2) + 3/2 Pi 
This result. was obtained by first computing the indefinite integral 
> int( 1/(5+3#cos(x)), x ); 


1/2 arctan(1/2 tan(1/2 x)) 
then by taking limits, given the following information about the discontinuities 


> readlib(discont) (",x); 
{2 Pi _Z1 + Pi} 


A second example showing that Maple can prove that this integral diverges, instead of 
returning unevaluated — meaning Maple cannot do it 


> int(1/x*2,x=-1..1); 

infinity 
Also, the addition of the CauchyPrincipalValue option to int handles 
> int( 1/x, x=-1..1, CauchyPrincipalValue) ; 


0 
> int( exp(-t)/t, t=-2..infinity, CauchyPrincipalValue ); 


- Ei(2) 


8.3. New and Enhanced System Facilities 


8.3.1 Type declarations for parameters 


Procedures accept type declarations as illustrated by the following example. The purpose of 
this facility is to encourage better type checking by making it easy to write and efficient. 


proc( a:t ) 


end 
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is equivalent to writing 


proc(a) 
if nargs > 0 and not type(a,’t’) then ERROR(message) fi; 


end 


Note, the type checking is not static, it is done when the procedure is called. The error 
message generated automatically is illustrated by this example 
> f := proc(n:integer) if n < 0 then 0 else n fi end: 


> f(x); 
Error, f expects its 1st argument, n, to be of type integer, but received x 


8.3.2 D extended to computation sequences and programs 


The D operator can now compute partial derivatives of functions which are defined as proce- 
dures. This is known as automatic differentiation. Consider 


> £ := proc(x) local t1,t2; 
> CL = 272: 

t2 := sin(x); 
3¥#t1*t2+24x¥*t 1-x¥*t2 


# The following computes the derivative of f wrt x, the first argument 


> 
> 
> end: 
> 
> Di] (£); 


proc(x) local tix,t2x,t1,t2; 
tix := 2x; 
ti := x72; 
t2x := cos(x); 
t2 := sin(x); 
3¥t1x*t2+3%t 1 et 2x+2¥t1xex+2*t1-t2x*x-t2 
end 
The reader can check that this really does compute the derivative of f by verifying 
that diff(f(x),x) - D(f)(x) = 0. The advantage of automatic differentiation is twofold. 
Firstly, it is much more efficient, in general, to represent a function as a program instead of 
as a formula. Secondly, it is more general, as functions can have conditional statements and 
loops. For example, given the array of the coefficients of a polynomial 6, we can represent 


the polynomial as a program in Horner form as follows. 


> £ := proc(x,b,n) local i,s; 

> s := 0; 

> for i from n by -1 to 0 do s := s¥*x+b[i] od; 
> s 

> end: 

> £(x,b,4); 


(((b[4] x + b[3]) x + b[2]) x + bli]) x + bf] 


> fx := D1] (£); 
fx := proc(x,b,n) 
local sx,i,s; 


sx := 0; 

s := 0; 

for i from n by -1 to 0 do sx := sx*xts; s := s*xt+b[i] od; 
sx 


end = 


8.3.3 Program optimization — optimize extended to procedures 


The optimize routine has been extended to optimize simple Maple procedures. Currently 
it does common subexpression optimization on computation sequences, i.e. procedures with 
assignment statements to local variables only, for example 
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> £ := proc(x) local t; t := x°2; 3*t*sin(x)+2*x*t-x¥*sin(x) end: 
> readlib(optimize) (f); 
proc(x) local t,t1; t := x72; ti := sin(x); 3*t*t1i+2*x*t-x*t1 end 


8.3.4 Multiple libraries — readlib takes multiple pathnames 


Maple automatically reads code from the Maple library using the readlib function. The 
global variable libname specifies the location of the Maple library. This can now be assigned 
a sequence of names. For example, if the the user assigns 
> libname := diri, libname, dir3; 
then when readlib(f) is executed, Maple will search for f first under dirl, and if un- 
successful, it will search under the Maple library, and if unsuccessful under dir3. This means 
users can have their own development library and have it override the Maple library. 


8.3.5 Assume facility 


One of the deficiencies of Maple and other systems is in handling problems which contain 
symbolic parameters where the result depends on the domain or range of values that the 
symbolic parameter(s) takes. For example, consider 


> Int( exp(-a*t)*ln(t), t=0..infinity ); 
infinity 
| 
| exp(- a t) 1n(t) dt 
| 
/ 
0 
The answer to this integral depends on the value of the parameter a. If a is real and 
positive, the answer is finite. If a is real and non-positive, the answer is infinite. How can 
the user tell Maple that a is real and positive? The solution adopted in Release 2 is to insist 
that the user state the assumption about a, i.e. that a is real and positive, as follows 
> assume(a>0); 
The result of this assumption is that the variable a has been assigned a new value which 
prints as a~ which Maple knows is real (implicitly) and positive. 


>a; 


a 


> about(a); 
Originally a, renamed a~: 
is assumed to be: RealRange(Open(0) , infinity) 


The assume facility is presently being integrated into Maple. In particular, Maple can now 
compute the definite integral above because Maple can determine that a is positive because 
signum(a) returns 1. 

> signum(a) ; 
1 


> int( exp(-a#t)*ln(t), t=0..infinity ); 


> is(at1>0); 
true 
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A deficiency of the assume facility is that the user needs to know what to assume in 
order to get answers out of facilities like definite integration. We are presently looking at 
alternative solutions including prompting the user for help or returning a solution along with 
the assumptions made. 


8.3.6 Automatic complex numeric arithmetic 


Complex arithmetic in Maple V had to be done with the evalc function. Now, complex 
numeric arithmetic is automatic, for example 


> (2+1/3)*2; 
35/9 + 4/3 I 


And complex floating point arithmetic 
> sin( (2.0+1/3)72 ); 


- 1.378776230 - 1.294704958 I 
The evalc function is now used only for exact symbolic expansions e.g. 


> evalc(exp(2+Pi/3*I)) ; 
1/2 
1/2 exp(2) + 1/2 I exp(2) 3 


> evalc( sqrt(atI*b) ); 


2 2 1/2 1/2 2 2 1/2 1/2 
(1/2 (a +b) + 1/2 a) + I csgn(b - I a) (1/2 (a +b) - 1/2 a) 
> assume(a>0); 
> assume(b>0) ; 
> evalc(sqrt(atb*I)); 
2 2 1/2 1/2 2 2 1/2 1/2 
(1/2 (a +b” ) + 1/2 a7) + I (1/2 (a~ +b” ) - 1/2 a7) 


8.3.7 Arrow operators now accept if statements 


The arrow operators, e.g. x -> x/(x°2-1); have been extended to handle if statements so 
that they can be used to define piecewise functions. E.g. x -> if x < 0 then 0 elif x < 
1 then x else 0 fi; 


8.3.8 I/O facilities and changes 


e printf: output a given expression with format (as in the C printf, with the added 
option of %a for algebraic), and also includes the C escapes: \b for backspace, \t for 
tab, \n for newline, e.g. 


“d\tvar = %s\tresult = %a\n‘, 10, y, 3*y*2+1); 
y result = 3*y"2+1 


> printf(‘n 
n= 10 var 


e sscanf: decode a string according to a format (inverse of printf, identical to the C 
sscanf function). 


input (tty and read) and output (1print and save) of numbers in scientific E/e nota- 
tion e.g. 1.2e+3 versus Float(12,2) or 1.2*10°3 


readline: read a line of an arbitrary file and return it as a string 
® parse: parse a string into a Maple expression 


e readdata: reads a file of numerical data arranged in columns 
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e output (lprint and save) of powers using ~ instead of ** 


8.3.9 Restart facility 


The restart statement will clear all user and system variables and reinitialize the kernel, 
essentially allowing you to start your Maple session from scratch without having to exit 
Maple. 


8.3.10 errorbreak 


Maple has added some user control for how Maple responds on encountering errors when 
reading in a file. This is provided by the interface option errorbreak. 


e With interface(errorbreak=0); the entire input file is read and all errors that are 
encountered are displayed. This was the Maple V behaviour. 


e With interface(errorbreak=1); processing stops after the first syntax error, thus 
avoiding nonsensical multiple syntax error messages that are often consequences of the 
first one. This is the new default behaviour. 


e With interface(errorbreak=2); processing stops after the first error, whatever it is, 
i.e. syntax (parse time) or semantic (execution) errors. 


8.4 Algorithmic and Efficiency Improvements 


8.4.1 Integration of rational functions 


The Trager-Rothstein algorithm computes the integral of a rational function f(x) over a field 
K as a sum of logs over K(a), an algebraic extension of K of minimial degree. E.g. consider 
the following rational function in Q(z) where note the denominator is irreducible over Q. 


625+ 823462? —42-—4 
26 +273—277+4+1 





f(z) = 
Applying the Trager-Rothstein algorithm yields the following result 


S> (l+a)In(z? - az +1) 


a?—2=0 


That is, there are two logs corresponding to two roots +2 of the polynomial a? — 2 = 0. 
Hence, Maple gets the following nice result 


> int(f,x); 


1/2 1/2 1/2 3 1/2 
(1 +2 ) In(x - 2 x+1)+ (i -2 +) Inte +2 x + 1) 

The first improvement in Release 2 is in handling coefficient fields which are function 
fields, and not just number fields. The input rational function could involve parameters as 
well as numbers. This is because of the extension of evala to handle algebraic functions and 
not just algebraic numbers. A second improvement due to Lazard and Trager allows for more 
efficient computation of the terms inside the logarithms. A third improvement by Rioboo 
improves the presentation of the result. An example of the latter is 


3 
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> £ := (x74-3%x72+6) / (x76-5#x744+54x7244) ; 


> int(f,x); 
The Trager-Rothstein algorithm (Maple V) obtains this result 
4/2 IT inte +I1x +3x-+21) - 1/21 inte ~Ix -3x+2 1) 
Rioboo’s algorithm (Release 2) yields this real result which has no new poles 


arctan(- 3/2 x + 1/2 x + 1/2 x) + arctan(x ) + arctan(x) 


8.4.2 Bivariate polynomial GCDs over Z 


We have implemented the dense modular method of Collins for computing GCDs of polyno- 
mials with in two variables with integer coefficients. This improves the efficiency of simpli- 
fying rational expressions in two variables. The method works by computing many GCDs 
of polynomials in one variable modulo different primes and combining these solutions using 
polynomial interpolation and the Chinese remainder theorem. Note, in Maple V, we im- 
plemented this method for univariate polynomials. Also, this is a dense method; for sparse 
polynomials in many variables it is not a good method. Comparing this with the default EEZ- 
GCD algorithm in Maple V, we obtain the following improvement. We created three dense 
bivariate polynomials g,a,b each of degree n with random n digit coefficients and computed 
the gcd(a*g, b*g). 





n 10 15 20 
Maple V_ 65.0 1,877 18,763 
Release 2 106 33.7 393.7 





8.4.3 Factorization over algebraic number fields 


We have implemented Lenstra’s algorithm for factorization of polynomials over an algebraic 
number field. For some polynomials it runs considerably faster than the default algorithm. 
But it depends very much on the polynomial. 


8.4.4 Faster numerical evaluation of the arctrig functions 


The running time of the algorithms used for the numerical evaluation of e” and In(z) at high 
precision is O(n?loggn) for n digits. This is O(loggn) integer multiplications each of which is 
O(n”). The algorithms for the trig functions sin, cos, and tan use repeated argument reduction 
to obtain an O(n5/2) algorithm. The same idea is now used by the arctrig functions, e.g. 


8.4.5 Optimized the evaluation of numerical arithmetic 


Evaluation of numerical expressions like n-i+1, 2*n-2*j+2, has been sped up to avoid the 
creation of an intermediate data structure. This is a significant improvement in the cases 
of small integers, e.g. array subscript calculations. Consider the following two routines for 
computing the Fibonacci numbers F'(n) 


> Fi := proc(n) if n < 2 then n else Fi(n-1)+Fi(n-2) fi end; 
> F2 proc(n) F2(n-1)+F2(n-2) end; F2(0) := 0; F2(1) := 1; 


46 Nieuwsbrief 








Digits sin! 7/4 tan! 7/4 
Maple V_ Release 2 | Maple V_ Release 2 | 
500 25 3.0 23 2.5 
1000 176 11.0 148 11.2 
2000 1245 50.4 1038 45.2 














The improvement in the arithmetic that occurs in F(n — 1) + F(n — 2) yields 





Maple V Release 2 
n Fl F2 Fl F2 


16; 1.05 0.73! 0.78 0.47 
18 | 2.738 217} 217 1,32 
20 | 7.30 5.20| 5.92 3.55 
22 | 20.17 13.45 | 15.20 8.53 

















8.4.6 Floating point solutions to polynomials 


We have rewritten the code for computing floating point approximations to roots of polyno- 
mials. The new code will increase the intermediate precision so that the roots found are fully 
accurate, i.e. accurate to 0.6 ulps. We have also implemented a routine in the kernel in C 
to use hardware floating point arithmetic. If the roots found using this routine are accurate 
then the fsolve routine will execute faster, as illustrated in the following examples. T(n, x) 
is the nth Chebyshev polynomial of the first kind which has all real roots. F(n,zx) is the 
nth Fibonacci polynomial which has all complex roots except 0. The jump in times indicates 
where hardware precision is no longer sufficient. 





Maple V Release 2 

nD) be) FC) | 2) Fae) 
10} 168 1.86; 045 0.23 
20; 8.95 11.43] 2.17 1.43 
30 | FAIL 36.72] 6.08 3.47 
40 FAIL | 34.22 = 7.57 
50 106.2 50.45 

















8.4.7 Character tables for Sn 


In computing the character table for S, the symmetric group on n letters we are now making 
use of the conjugate partitions and symmetry in the table. The improvement is better than a 
factor of 2 for large n. Note the table has dimension p(n) by p(n) where p(n) is the number of 
partitions of n so requires exponential space in n to store. The data here are for computing 
combinat [character] (n) 


Secs: 3: 
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4T 
n | Maple V | Release 2 
8 6.6 3.5 
10 26.0 12.6 
12 102.0 43.4 
14 446.8 160.3 








New Library Functions 


AFactor, AFactors: absolute factorization over the complex numbers 
argument: complex argument of a complex number 

assume: for making assumptions about variables 

ceil: ceiling function 

csgn: complex sign 

discont: compute a set of possible discontinuities of an expression 
floor: floor function 

fourier, invfourier: symbolic fourier transform and inverse 

ilog, ilogi0: IEEE integer logarithm function 

Indep: test algebraic extensions over Q for independence 

invfunc: table of inverse functions 

is, isgiven: for testing properties against the assume database 
makehelp: for making a Maple help page from a file of text 

maxorder: an integral basis for an algebraic number or function field 
Norm: computes the norm of the algebraic number or function 
Primfield: computes a primitive element over a given algebraic number field 
powmod: computes a(x)" mod 6(z) using binary powering 

ratrecon, Ratrecon: rational function reconstruction (Euclidean algorithm) 
readdata: reads a file of numerical data arranged in columns 

spline: computes a natural cubic spline from a set of data 

split: splits a polynomial into linear factors over its splitting field 
sqrfree: a square free factorization of a multivariate polynomial 
symmdiff: symmetric difference of sets 

tutorial: on-line tutorial facility 

Trace: computes the trace of the algebraic number or function 


unload: unload a Maple library function 
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8.6 Library Packages 


The package mechanism has been extended to support an initialization routine. The initial- 
ization routine should be placed either in the package as the routine package[init], or as 
the routine ‘package/init‘ in the library in the file ‘package/init.m‘; This routine is 
executed by with(package) ; It is up to the initialization routine to decide what to do if it 
is called a second time. Below is a list of the new packages and modifications to existing 
packages. 
8.6.1 numapprox (Numerical Approximation Package) 

e pade: computes a Pade rational function approximation 

e infnorm: compute the (weighted) infinity norm of an approximation 


e minimax: compute a (weighted) minimax numerical approximation 


e remez: Remez algorithm for minimax rational function approximation 


chebpade: compute a Chebyshev-Pade rational function approximation 


chebyshev: compute a Chebyhev series numerical approximation 
e confracform: convert a rational function to a continued fraction 
e hornerform: convert a polynomial or rational function to Horner form 


The call minimax (f(x) ,x=a..b, [m,n] ,w) uses the Remez algorithm to compute the best 
minimax rational function with numerator of degree < m and denominator of degree < n 
(Note: n can be zero, meaning a polynomial) to the function f(x) on the range (a, 6]. The 
optional 4th paramter w specifies a weight function (default 1). For example, here is a 
[4,3] rational approximation to erf(x) on [0,3] which we’ve rounded to 7 decimal digits. The 
approximation is accurate to 4 digits on [0,3]. 

> evalf (minimax (erf (x) ,x=0..3,[4,3]), 7); 


- .000088319 + (1.132571 + ( - .4507723 + (.1724417 - .0209750 x) x) x) x 


1. + ( - .3706125 + (.4051100 - .04578747 x) x) x 


8.6.2 DEtools (Differential Equation Tools Package) 


A package for plotting differential equations including systems of ODE’s and some special 
PDE’s, allowing the associated vector field and solution curves for different initial conditions, 
to be plotted. 


e DEplot: plot a single ODE or a system of first order ODEs 
e DEploti: plot a first order ODE 
e DEplot2: plot a system of two first order ODEs 


e PDEplot: plot the surface of a first order, quasi-linear PDE of the form 
P(x,y,u) * D{1](u) (x,y) + Q(x,y,u) * D[2](u) (x,y) = R(x,y u) 


e dfieldplot: plot a field plot for a single or two first order linear ODEs 


e phaseportrait: plot a phase plot for a single or two first order linear ODEs 
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8.6.3. Gauss 


A package which allows one to create ”domains” of computation in a similar way that the 
AXIOM system allows. The advantage is that routines for computing with polynomials, 
matrices, and series etc. can be parameterized by any coefficient ring. For example, given a 
ring R, the Gauss statement 
> P := DenseUnivariatePolynomial(R,x): 

creates a dense univariate polynomial domain P which contains all the routines needed to 
do arithmetic with polynomials over R. The advantage is that the the code supplied by the 
DenseUnivariatePolynomial function works for any coefficient domain R which is a ring, i.e. 
has operations +, -, *, etc. It is expected that many Maple library routines will be recoded 
in Gauss over the next few versions. The initial implementation includes basic facilities for 
computing with polynomials, matrices and series over number rings, finite fields, polynomial 
rings, and matrix rings. 


8.6.4 GaussInt (Gaussian Integer Package) 


This is a set of routines for computing with Gaussian (complex) integers, i.e. numbers of the 
form a + b I where a and b are integers. Included are routines for computing the GCD of 
two Gaussian integers, the factorization of a Gaussian integer, and a test for primality, e.g. 


> GIfactor(13); 
(3 +21) (3 - 2 I) 


8.6.5  linalg (Linear Algebra Package) 


e blockmatrix: utility routine to create a block matrix 


e gausselim, gaussjord, rowspace, colspace: now handle matrices of floats, com- 
plex floats, and complex rationals 


eigenvals, eigenvects: now handle matrices of floats and complex floats 


entermatrix: utility routine for entering the entries of a matrix 
e randvector: utility routine to create a random vector 
e ratform: (synonym for frobenius) computes the rational canonical form 


e Wronskian: computes the Wronskian matrix 


8.6.6 numtheory: (Number Theory Package) 


e cfrac: modified to handle different forms of continued fractions real numbers, series, 
polynomials. 


e cfracpol: computes continued fractions for the real roots of a polynomial 
e kronecker: Diophantine approximation in the inhomogeneous case 

e minkowski: Diophantine approximation in the homogeneous case 

e nthnumer: returns the denominator of the nth convergent 

e nthdenom: returns the numerator of the nth convergent 


e nthconver: returns the nth convergent of a simple continued fraction 
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e thue: compute the resolutions of Thue equations or Thue inequalities 


e sq2factor: factorization in the UFD 2Z[V2] 


8.6.7 padic (P-adic Number Package) 


A package for computing P-adic approximations to real numbers, e.g. 


> with(padic): 
> evalp(sqrt(2),7); 


2 3 4 5 6 T 8 9 
3+7+27 +67 +7 +27 +7 +27 +47 «++0(7) 


> evalp(exp(3) ,3); 


2 3 4 6 8 9 
1+3+3 +23 +23 +3 +3 +03) 


8.6.8 plots (Plots Package) 


e animate: animate one or more functions in 2 dimensions, e.g. 
a) sine - cosine waves : 
animate({sin(x*t) ,cos(x#t)},x=-2*Pi..Pi,t=0..2*Pi); 
b) a polar coordinate unraveling 


animate ([sin(x*t) ,x,x=-4..4],t=1..4,coords=polar,numpoints=100,frames=100) 


e animate3d: animate one or more curves or surfaces in 3 dimensions, e.g. 


animate3d(cos(t*x)*sin(t*y), x=-Pi..Pi, y=-Pi..Pi,t=1..2); 


conformal: routine is typically 3 times faster 


contourplot: contour plot of a 3 dimensional surface, e.g. 
contourplot(sin(x*y) ,x=-Pi..Pi,y=-Pi..Pi,grid=[30,30]); 

e densityplot: density plot, e.g. 

densityplot(sin(x*y) ,x=-Pi..Pi,y=-Pi..Pi,grid=[30,30]); 


e fieldplot: plot of a 2 dimensional vector field, e.g. 
vfield := [x/(x"2+y°2+4)*(1/2) ,-y/(x*2+y°2+4)*(1/2)]: 
fieldplot (vfield,x=-2..2,y=-2..2); 


e fieldplot3d: plot of a 3 dimensional vector field, e.g. 
fieldplot3d([2*x,2*y,1] ,x=-1..1,y=-1..1,z=-1..1, 
grid=[10,10,5] ,arrows=THICK,style=patch) ; 


gradplot: plot the gradient vector field of a 2 dimensional function 


gradplot3d: plot the gradient vector field of a 3 dimensional function, e.g. 
gradplot3d( (x"2+y°2+z°2+1) (1/2) ,x=-2..2, y=-2..2,z=-2..2); 

e implicitplot: plot a curve defined implicitly by F(z, y) = 0, e.g. 

implicitplot( x°2 + y°2 = 1, x=-1..1,y=-1..1 ); 
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implicitplot3d: plot a surface defined implicitly by F(z, y,z) = 0, e.g. 
implicitplot3d( x*2ty72+z°2 = 1, x=-1..1, y=-1..1, z=-1..1 ); 


logplot: log plot, e.g. logplot (107 (20*exp(-x)) ,x=1..10); 


loglogplot: log log plot, e.g. loglogplot (10*(20*exp(-x)) ,x=1..100); 


odeplot: plot the solution of an ODE from dsolve 
p := dsolve({D(y)(x) = y(x), y(0)=1}, y(x), numeric): 
odeplot(p, [x,y(x)],-1..1); 


e polygonplot: draw one or more polygons in 2 dimensions 
ngon := n -> [seq([ cos(2*Pi*i/n), sin(2*Pi*i/n) ], i = 1..n)]: 


polygonplot (ngon(8)); 
e polygonplot3d: draw one or more polygons in 3 dimensions 


e polyhedraplot: for plotting tetrahedrons, octahedrons, hexahedrons, dodecahedrons, 
icosahedrons e.g. 


polyhedraplot([0,0,0] ,polytype=dodecahedron) ; 


surfdata: plot a surface input as a grid of data points 
cosdata := [seq([ seq([i,j,evalf(cos((itj)/5))], i=-5..5)], j=-5..5)]: 
F y= (xsy) -> x2 + 72: 


surfdata( cosdata, axes=frame, color=F ); 


textplot: for plotting labelled points in 2 dimensions, e.g. 
textplot([1,2,‘one point in 2d‘],align={ABOVE,RIGHT}) ; 


e textplot3d: for plotting labelled points in 3 dimensions, e.g. 
textplot3d({[1,2,3,‘one point in 3d‘],(3,2,1,‘second point in 3d‘]}); 


8.6.9 networks: (Graph Networks Package) 


A package of routines for creating, manipulating and computing properties of graphs, i.e. 
graphs with edges and vertices. Primitives exist for computing properties such as edge con- 
nectivity, biconnected components, and edge or vertex set induced subgraphs. The graphs 
may have loops or multiple edges, both directed and undirected, and arbitrary weights can be 
associated with both edges and vertices. Graphs can be drawn as a picture (a plot) for display. 
For example, here we create the Petersen graph and compute and factor its characteristic 
polynomial. 


> with(networks) : 
> factor(charpoly(petersen() ,x)); 


4 5 
> 3) G@+2). G@- 
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The spanpoly(G,p) function computes a polynomial for the undirected graph G where 
p is a probability of failure associated with each edge. When G is connected, this is the 
all-terminal reliability polynomial of G. For example, the all-terminal reliability polynomial 
for the complete graph on 4 vertices 


> spanpoly(complete(4) ,p); 


3 2 3 
-~p (- 16+ 33 p- 24p +6p) 


8.7 Share Library 


8.7.1 Contributions from Maple users 


The share library is a repository of applications codes, contributed primarily by Maple users. 
It is distributed with each version of Maple in both src and Maple ”.m” format. The share 
library is being reorganized. Beginning with this version, it also includes 


e TeX/LaTeX introductory documents on using Maple 

e TeX/LaTeX support documentation for applications routines 

e updates to library routines, i.e. bug fixes, efficiency improvements etc. 
e examples of applications of Maple to problem solving 


We have also endeavoured to ensure that each contribution comes with on-line Maple style 
help files, and Maple style test files to help us ensure that the code is working correctly. 


8.7.2 How to access code from the share library from inside Maple 


Users can access routines in the share library by reading them in directly from the files using 
the read command. Also, where appropriate, one can now access them directly from inside 
Maple using the with command. 

The command with(share) ; looks to see if the share library exists, i.e. it is adjacent to 
the Maple library (it simply tries to read a file from it). If successful, it assigns the global 
variable sharename to the path where it found the share library, and then it executes the 
statement 


> libname := libname, sharename; 


Now Maple’s readlib command and with command will automatically search the share 
library after first searching the Maple library. If unsuccessful, the user can locate the share 
library and set the libname variable directly. For example, the fit routine in the share library 
can be used to fit a formula to some data in the least squares sense. 

> with(share): # locate the share library 


> readlib(fit): # now finds the fit routine from the share library 
> £it([-1.5,0.5,1.5],[0,-2,-1.5], axt*2+bet+c, t); 


2 
-4999999997 t - .5000000001 t - 1.874999999 


Documentation for the code is in the file, ie. after reading the file, ?7fit will print out a 
regular Maple style help file. Note, in many cases, additional documentation is included in 
TeX/LaTeX files. See ?share,contents for a summary of the contents of the share library. 


“~) 
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8.7.3. Electronic distribution version of the Maple share library 


The share library will also be distributed electronically, firstly, over the internet by using 
anonymous ftp, and secondly, by electronic mail, with an electronic mail server. The elec- 
tronic version is be updated periodically, about once every 3 months, with new contributions, 
corrections to previous contributions etc. Unlike the Maple version, it contains a directory of 
updates to the Maple library, primarily simple bug fixes. For directions on how to get code 
using anonymous ftp, electronic mail, and how to install it and use it, see ?share for details. 
8.7.4 Policy of the contributions to the share library 


The code in the share library has been contributed by the authors to be made freely available 
to all Maple users. It is not to be sold for profit or changed without the author’s permission. 
Authors who contribute code to the share library will be asked to write Maple style on- 
line help documentation and a test file. Authors will also be asked to sign a non-exclusive 
agreement to let us distribute your code, and to change it to maintain it. 


8.8 Miscellaneous 


e the catenation operator . now expands in ‘‘.(a,b,c).(1..3) to yield al, a2, a3, bl, 
b2, b3, cl, c2, c3 


e expand of @ and @@ now handled e.g. expand((sin@@)(2)(x) ); ==> sin(sin(x)) 
e evalf of @ and @@ now handled e.g. evalf((sin@@2)(2)); ==> .7890723436 

e anames(t) returns all names of type t 

e ali] := x; is dissallowed if a is assigned an object other than an array or table 

e £()[1] is now allowed; it returns unevaluated 

e Changes to fortran and C 


— use appendto instead of writeto to append output to a file instead of overwriting 
it 


both now translate subscripts e.g. a[i-2*j] where the 2 here is understood to be 
an integer 


both now convert symbolic constants e.g. Pi to floating point 


~ optional parameter digits = n for specifying the precision to be used in convert- 
ing constants to floating point 


— C now also breaks up large expressions which can cause some compilers to break 


Changes to dsolve/numeric 


The output of dsolve({...},{...},numeric) is a procedure f where f(z) returns 
a tuple of values. E.g dsolve({...},{y(x),z(x)},numeric); returns f a procedure 
where f(1) returns x = 1, y(x) = y(1),2(x) = z(1) The odeplot routine in the plots 
package can be used to plot the result. 


mellin speeded up by a factor of 4-30 


e ifactor has been sped up for the case of a perfect power 
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series (f(x) ,x=a); now allows tooa 
> series(1n(1-x) ,x=-infinity) ; 


1 1 1 1 1 

~ In(- 1/x) - 1/x - ---- - ---- - ---- - ---- + 0(----) 
2 3 4 5 6 
2x 3x 4x 5 x x 


signum(f), where f is a symbolic constant will use interval arithmetic at Digits pre- 
cision to determine the sign of f, e.g. signum(Pi-4*sin(1)); 


the for loop clauses by, to, from and while can now appear in any order 
Changes to automatic simplifications 
— the 1n(1/x) ==> ~1n(x) automatic simplification has been turned off. In(a*n) ==> n J 
is now done iff signum(a) = 1, ie. x j 0 is provable 


— the ln(exp(x)) ==> x automatic simplification has been turned off. In(exp(x)) ==> x 
is now down iff Im(z) returns 0, i.e. x is provable real 


— the automatic simplification sqrt (x2) ==> x (incorrect for x < 0) is being fixed 
for the next version (this change is more difficult than the others) 


— the arcsin(sin(x)) ==> x arctrig simplifications have been turned off 
~ inverse functions simplifications have been added e.g. sin@@(-1) ==> arcsin 
— trig simplifications for multiples of 7/8, 7/10, 7/12 have been added 

> sin(Pi/8), sin(Pi/10), sin(Pi/12); 


1/2 1/2 1/2 1/2 1/2 
1/2 (2-2 .) , 1/45 - 1/4, 1/2 (2-3 ) 


given a file of text, the makehelp function makes a Maple TEXT object out of it, i.e. 
an on-line help page. I.e. you can make your own help pages which can be displayed 
using the ? command using this. 


Additional simplifications for Dirac and Heaviside, and the arc trigonomentric func- 
tions are provided through the simplify facility, e.g. 
> simplify( Heaviside(x-1)*Heaviside(x-3) ); 


Heaviside(x - 3) 
> simplify( F(x+3)*Dirac(x) ); 

F(3) Dirac(x) 
> simplify( arctan(x)+arctan(1/x) ); 


1/2 signum(x) Pi 


] 


or: UO 


ee ee ee ae a oe 


— PS | ee eS es ee eng are 


9 | Programmeren in Maple 


W.T. Goes 
Knowledge-Based Engineering Group 
Faculty of Technical Mathematics and Informatics 
Delft University of Technoligy 
The Netherlands 


9.1 Wisselvalligheid 


Stel we willen een routine in Maple schrijven die bij gegeven invoer 1/a, 2/a, 3/a, ..., de 
generalisatie k/a moet vinden. Bij het maken van deze routine kunnen we gebruik maken 
van de in Maple aanwezige functies ’whattype’, ’op’ en ’nop’. 

De functie ’whattype’ geeft de volgende antwoorden: 


e whattype(1/2); levert in Maple ’fraction’. 
e whattype(2/a); levert in Maple ’*’. 
e whattype(1/a); levert in Maple’ *’. 


De expressies 1/2, 2/a en 1/a worden in de wiskunde als breuken gezien. Maple echter 
maakt er drie verschillende types van (breuk, vermenigvuldiging, exponent). Ook de functies 
‘op’ en ’nops’ vertonen dit gedrag, omdat ze rechtstreeks in verband staan met de functie 
*whattype’. Deze wisselvallige classificaties komen niet overeen met de verwachtingen van 
een gebruiker. 

In dit artikel worden de volgende vragen beantwoord: 


e Wat is de oorzaak van deze wisselvalligheid? 
e Hoe kan dit gedrag door de producenten verbeterd worden? 


e Hoe kan dit gedrag door de gebruiker verbeterd worden? 


9.2 De oorzaak 


Computers werken doordat reeksen instructies de taak omschrijven die uitgevoerd dient te 
worden. Op het allerlaagste niveau van een computer bevinden zich eenvoudige en elementaire 
instructies. Deze instructies zijn gemakkelijk in hardware te implementeren, maar moeilijk 
door mensen te gebruiken. 

We kunnen vervolgens een machine verzinnen welke een gebruikersvriendelijkere instruc- 
tieset heeft en alle technische details van het vorige niveau verbergt. Deze machine is door 
de gebruiker beter te bevatten. Als nog een hogere abstractie gewenst is, kan het voorgaande 
proces nog een aantal malen toegepast worden. Afbeelding 1 laat deze denkwijze zien. 
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Deze techniek, waarbij gebruik gemaakt wordt van denkbeeldige machines, heeft een zeer 
ruim toepassingsgebied. Hierbij moet men denken aan communicatienetwerken en comput- 
ers, maar ook aan tv’s en wasmachines. De gebruikers van deze apparaten behoeven zich 
niet te bekommeren om details, maar kunnen gebruik maken van de beschikbare instructies 
(knoppen) op een gemakkelijk te bedienen niveau. 

De verschillende virtuele machines die bijvoorbeeld op een computer voor kunnen komen, 
zijn getoond in de figuur 2. 

Ook voor de gebruiker van Maple is het van wezenlijk belang om abstractieniveau’s te 
onderscheiden: de wiskundige functies (zoals diff en limit) en de interne representatie van 
Maple daartoe. 

Daarom zouden we bij de in Maple aanwezige programmeertaal ook dergelijke abstracties 
verwachten: een scheiding tussen de voor de computer efficiénte interne structuur en de voor 
de gebruiker makkelijk te manipuleren hogere orde taal. Helaas is Maple hier niet consistent 
in. Dit is getoond in figuur 3. 

Zoals in figuur 3 uitgebeeld is, ziet de gebruiker de taal Maple als een mengeling van de 
Internal representation’ en een "High Level Language’ (HLL). Het blok HLL schermt niet de 
gehele interne representatie van Maple voor de gebruiker af. 

Laten we nog even kijken naar het voorbeeld uit de inleiding. De resultaten die Maple 
geeft bij de instructie ’whattype’ zijn een weergave van de structuur die de computer zelf 
hanteert. Als we de functie ’whattype’ onafhankelijk van deze interne representatie zouden 
maken, kunnen we een functie maken die intuitief beter is. Bijvoorbeeld: 


e Whattype(1/2); zou moeten leveren ’fraction’ 
e Whattype(2/a); zou moeten leveren fraction’ 


e Whattype(1/a); zou moeten leveren ‘fraction’ 


e Whattype(a* b); zou moeten leveren ’* ’ 


( 
( 
( 
e Whattype((x+2)/(y+2)); zou moeten leveren ’fraction’ 
( 
( 


Whattype(a*c*g); zou moeten leveren ’*’ 


De kennis van Maple en van de wiskundige op het terrein van de wiskunde komen nu 
meer overeen. 

In het volgende hoofdstuk zal eerst een korte omschrijving gegeven worden van Maple’s 
interne representatie. Na het lezen hiervan kunt u de overeenkomst zien tussen deze rep- 
resentatie en bijvoorbeeld de functie ’whattype’. Hierna worden mogelijkheden gegeven om 
Maple’s wiskundige kennis bij te schaven. 


9.3. Interne representaties 


Alle in Maple aanwezige interne datastructuren hebben dezelfde algemene structuur. Deze 
stuctuur bestaat uit een reeks velden. Het eerste veld in deze reeks geeft het type van de 
structuur aan. De velden die daarop volgen, bevatten de bij de structuur behorende data. 
Het volgende voorbeeld toont deze algemene structuur: 


e TYPE — DATA — DATA — DATA — ENZ. 


De structuren ‘Product’ en Sum of terms’ zullen verder uitgewerkt worden. Een ver- 
menigvuldiging van de gebruiker is intern aanwezig d.m.v. een ’Product’-structuur. Deze 
representatie is als volgt: 


br eink oak 
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e PROD — EXPR — EXPO — EXPR — EXPO — ENZ. 


Het eerste veld geeft het type aan. Deze structuur representeert het product van een reeks 
subexpressie-exponent combinaties. De datavelden geven de expressies met hun exponenten 
weer. Een voorbeeld van een expressie representeerbaar in deze structuur is (a+ 1)*b* (c+ 
e+3)°d. Deze expressie wordt intern met behulp van de ’Product’-structuur op de volgende 
manier weergegeven: 


e PROD — (a+ 1) — b— (c+e+3) —d 
De tweede structuur Sum of several terms’ is als volgt: 
e SUM — EXPR — FACTOR — EXPR — FACTOR — ENZ. 


Deze structuur representeert de som van een reeks factor-subexpressie combinaties. Een 
expressies van dit type is 2* (a+ z)+3% (b+ <2) en deze wordt op de volgende manier intern 
gerepresenteerd: 


e SUM — (a+2z) —2—(b+2)—3 
De expressie 2/a wordt gerepresenteerd met een subexpressie 1/a en een factor 2: 
e SUM — 1/a — 2 


Bij voortzetting van deze methode zou 1/a gerepresenteerd moeten worden met een subex- 
pressie 1/a en een factor 1: 


e SUM — 1/a—1 


Daar dit tot een oneindige herhaling leidt, representeert Maple 1/a op een andere manier. 
De expressie 1/a wordt door Maple gerepresenteerd met behulp van de ’Product’-structuur 
met een subexpressie a en een exponent —1: 


e PROD — a— -1 


Hierdoor is de oneindige herhaling vermeden. Het hier aangegeven onderscheid tussen 
2/a en 1/a wordt ook gemaakt door de functie ’whattype’. Het wisselvallige gedrag van deze 
functie is dus een rechtstreeks gevolg van de interne representatie. 

Een ander representatieprobleem komt voor de gebruiker te voorschijn bij functies die 
lijsten als uitvoer leveren. Deze lijsten worden gerepresenteerd d.m.v. het type ‘expression 
sequence’ (EXPSEQ). Dit type wordt intern gebruikt om ’tables’ en ‘arrays’ op te slaan. Het 
is echter niet mogelijk om een lijst van het type EXPSEQ weer rechtstreeks te gebruiken als 
parameter van een nieuwe functie. 

Als we bijvoorbeeld het resultaat van de ’op’ operatie aan een variabele ’L’ toekennen, 
krijgen we: 


e L := a,b,c,2; 


De actie ’nops(L)’ genereert vervolgens een foutmelding, omdat deze actie gezien wordt 
als: 


e nops(a,b,c,2); 
en teveel parameters heeft. De correcte aanroep in het huidige Maple is: 
e nops([L}); 


De rechte haken zijn een hulpmiddel om het gebruik van de interne representatie recht 
te zetten. Deze rechte haken forceren een array-representatie. Immers, als alle resultaten 
die van het type ’EXPSEQ’ zijn, vervangen worden door het type ‘array’, hoeft u zich als 
gebruiker geen zorgen meer te maken om de interne weergave van Maple. 
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9.4 Mogelijke veranderingen 


De in Maple aanwezige wiskundige functies vormen de kracht waaraan het programma zijn 
reputatie ontleent. De veranderingen die hieronder voorgesteld worden, laten het wiskundige 
gedeelte daarom natuurlijk volledig intact. De aanwezige programmeerfaciliteiten staan cen- 
traal in de hier volgende besprekingen. 

Om Maple met een beter programmeerbaar en voor de gebruiker consistenter systeem uit 
te rusten, bekijken we een tweetal mogelijkheden. 


e Wiskundige faciliteiten van Maple onderbrengen in een bestaande programmeertaal. 


e Maple uitbreiden met meer consistentere programmeerfaciliteiten, zodat een ’hogere 
orde’ situatie ontstaat. Deze oplossing reeds is geschetst. 


9.4.1 Maple onderbrengen in een bestaande taal 


Het doel van deze mogelijkheid is de wiskundige kracht van Maple te combineren met de 
programmeerkracht van bijvoorbeeld C++. Dit combineren kan op de volgende manier be- 
werkstelligd worden. Plaats de functionaliteit van Maple in een aantal verschillen DLL’s 
(Dynamic Link Library). Figuur 4 toont de wijze waarop de DLL’s gekoppeld worden. 

Onder OS/2 en Windows worden de hiervoor genoemde DLL’s toegepast. 

De basis van de structuur wordt gevormd door de units ’Built-in functions’ en ’Built- 
in types’. Deze units bevatten de zogenaamde ’Maple kernel’. Dit gedeelte definieert de 
basisoperaties (procedures en functies), welke in de sourcecode van bibliotheken en applicaties 
gebruikt kunnen worden. In de unit ’Libraries’ bevinden zich de door Maple meegeleverde 
wiskundige bibliotheken. Deze bibliotheken worden gegeven in de sourcecode van de taal 
waarin Maple opgenomen is. Op deze manier kunnen bibliotheekfuncties door de gebruikers 
aan hun eigen wensen aangepast worden. Dit concept wordt ook in de huidige versie van 
Maple toegepast, alleen wordt er dan gebruikt van de eigen Maple code. 

Met de gegeven basisunits en de bibliotheekunits kunnen nu applicaties geschreven wor- 
den. De programmeur van applicaties heeft nu de beschikking over interfacefaciliteiten zoals 
XView of Windows, alle mogelijkheden van bijvoorbeeld C++ of Pascal en natuurlijk de 
wiskundige capaciteiten van Maple. Om compatibiliteit met de huidige versies van Maple 
te garanderen, dient een applicatie reeds aanwezig te zijn. Deze applicatie interpreteert de 
huidige Maple code. Hiermee kunnen alle in Maple code geschreven toepassingen gebruikt 
worden. 


9.4.2 Maple uitbreiden met consistentere functies 


Het onderbrengen van Maple in een programmeertaal is een aardige oplossing, maar het 
bewerkstelligen van deze verandering ligt niet in het bereik van de gebruiker. De tweede 
mogelijkheid, die nu besproken gaat worden, kan op beperkte schaal door gebruikers van 
Maple uitgevoerd worden. 

De interne representatie van expressies in Maple is direct zichtbaar bij het toepassen van 
diverse functies. Het doel is nu deze functies aan te passen of te vervangen, zodat deze 
interne representatie aan het zicht onttrokken wordt. Deze mogelijkheid zal in meer detail 
in de volgende paragraaf uitgewerkt worden. Als voorbeeld wordt de functie *whattype’ 
herschreven. 


Ses a, ve a ae ee SE: 


+” 
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9.5 Aangepaste functies 


9.5.1 Overzicht van mogelijke veranderingen 


Zoals in de inleiding al aangegeven is, dient de functie whattype’ veranderd te worden. Deze 
verandering wordt hierna gedetailleerd uitgewerkt. De verbeterde functies beginnen met een 
hoofdletter. 

Bij een verandering van de functie ’whattype’ is een verandering van een aantal andere 
functies ook noodzakelijk geworden. De bestaande functies ’op’ en ’nops’ zijn afhankelijk 
van de resultaten van de huidge ’whattype’. Zo ook de nieuwe ’Op’, ’Nops’ en ’Whattype’. 
Bijvoorbeeld bij de expressie (c * b)/a dienen de resultaten te zijn: 


e Whattype((c*b)/a); zou moeten leveren ’fraction’ 
e Op((c*b)/a); zou moeten leveren ’array([c*b,a])’ 


e Nops((c*b)/a); zou moeten leveren ’2’ 


De veranderde operatie ’Op’ levert een array als resultaat. Het type EXPSEQ is hiermee 
volledig aan het zicht van de gebruiker onttrokken. Andere bestaande functies die ook op 
typen betrekking hebben en daarom ook voor verbetering in aanmerking komen, zijn: ’type’, 
*hastype’ en ’convert’. 

Een bestaande functie welke ook in aanmerking komt voor verandering, is de functie ’sort’. 
Bij het sorteren van expressies wordt in Maple geen gebruik gemaakt van een wiskundige 
ordening. Voorbeelden van de huidige situatie zijn: 


e sort([x+2,x+3]); levert in Maple ’[x+2,x+3] 
e sort([x+1,x+2]); levert in Maple ’[x+2,x+1]’ 


Omdat dit tegenstrijdig is met de verwachtingen, is ook een aangepaste functie ’Sort’ 
gewenst. 


9.5.2 whattype 


De functie ’whattype’ bepaalt het basistype van een expressie. Het gedrag van deze functie 
dient veranderd te worden zoals in de inleiding aangegeven is. De aanroepen van de nieuwe 
functie 'Whattype’ zijn met een hoofdletter geschreven, omdat het om een wezenlijk andere 
functie gaat. 


e Whattype(1/2); zou moeten leveren fraction’ 

e Whattype(2/a); zou moeten leveren fraction’ 

e Whattype(1/a); zou moeten leveren fraction’ 

e Whattype((x+2)/(y+2)); zou moeten leveren ‘fraction’ 
e Whattype(a* b); zou moeten leveren ’*’ 


e Whattype(a*c*g); zou moeten leveren ’*’ 


De functie nieuwe ’Whattype’ wordt hieronder verder uitgelegd. De programmacode 
wordt daarna gegeven. 

Bij het definiéren van de nieuwe functie worden aan de fend van interne definities van 
Maple drie gevallen onderscheiden. 
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e ’Product’-structuur. Deze structuur representeert het product van een reeks subexpressie- 
exponent combinaties. 


e ’Sum of several terms’-structuur. Deze structuur representeert de som van een reeks 
factor-subexpressie combinaties. 


e Alle resterende structuren. 


Als bij de interne representatie van een expressie geen gebruik gemaakt wordt van een 
van de eerste twee situaties dan wordt de originele ’whattype’ functie aangeroepen. In de 
andere twee gevallen worden de afwijkende situaties afgehandeld. 

Bij de ’product’-structuur is er sprake van het type ’fraction’ als er minstens een element 
voorkomt met een negatieve quotient. In alle andere gevallen wordt er gebruik gemaakt van 
de originele ’whattype’ functie. Op deze manier worden alle als breuk uitziende elementen ook 
als zodanig geclassificeerd. De routine voor dit gedeelte is gegeven onder de naam ’Quotient’. 

Bij de ’sum of several terms’-situatie ligt het iets ingewikkelder. Komen meerdere termen 
voor dan wordt het een expressie van het type ’+’. In het geval van een enkele term worden 
echter meerdere situaties onderscheiden: 


e Het expressie gedeelte bevat een negatieve quotient. 
e De factor is een breuk. 


e Alle resterende situaties 


In de twee bovenste gevallen wordt het type geclassificeerd als een ‘fraction’. In de nog 
resterende gevallen wordt het type bepaald door de originele functie ’whattype’. De routine 
voor de gehele vervangende functie draagt de naam ’Whattype’. Deze functie maakt gebruik 
van de functie Quotient’. 


9.5.3 Programmacode ’Quotient’ 


Quotient := proc() 

local _w_e, _wlist, -w_cnt, -.w_loop; 

_w_e := args(1]; 

-w list := disassemble(addressof(-w_e)); 

_w_cnt := 0; 

for _wloop from 3 by 2 to nops([-w_list]) do 
if pointto(-w_list[.w loop]) < 0 then _w_cnt := _w_cnt + 1 fi 

od; 

if w_cnt <>0 then RETURN(fraction) else 
RETURN(whattype(_w_e)) 

fi 

end 


9.5.4 Programmacode ’Whattype’ 


Whattype := proc() 
local _w_e, _w_list, w_type; 
_w_e := args[1]; 
-w list := disassemble(addressof(_w_e)); 
_w_type := _w_list/[1]; 
if w_type = 11 then 


ee ee Ee ee ee ee ee ee ee ee eee 


ee Pe ee ee ee ee ee oe 
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if nops([-w-list]) = 3 then 
if disassemble(_w_list[3])[1] = 3 then RETURN (fraction) 
elif disassemble(_w_list[2])[1] = 9 then 
if Quotient(pointto(_w_list[2])) = fraction then 
RETURN (fraction) 
ij 
RETURN(whattype(_w_e)) 
else RETURN(whattype(_w_e)) 
fi 
else RETURN(‘+*‘) 
fi 
elif .w_type = 9 then RETURN(Quotient(_w-e)) 
else RETURN(whattype(_w-e)) 
fi 
end 


9.6 Conclusie 


Wisselvallig gedrag van Maple is het gevolg van de interne weergave van expressies. De 
gebruiker wordt hierdoor afgeleid van zijn of haar wiskundige taak en dient zich te verdiepen 
in de door de computer gebruikte interne structuren. 

Het onzichtbaar maken van deze interne representaties zorgt ervoor dat de wiskundige 
inzichten van Maple meer over eenkomen met die van de wiskundige gebruikers. Frustraties 
die het gevolg zijn van onverwachte resultaten worden hiermee sterk verminderd. 

Dit aanbrengen van een hoger abstractie niveau is dus van wezenlijk belang en kan op 
kleine schaal door de gebruiker zelf uitgevoerd worden. 


9.7 Referenties 


Met dank aan prof. dr. H. Koppelaar voor de gegeven ondersteuning en voor het doorlezen 
van eerdere versies van dit artikel. 

B.W. Char, K.O. Geddes, G.H. Gonnet, e.a., Maple V Language Reference Manual, 
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Figuur 2: Concreet voorbeeld 
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Figuur 3: Situatie bij Maple 





Figuur 4: DLL structuur 


64 


Nieuwsbrief 


10 | Report on the Spring School on the 
Didactics of Computer Algebra! 


10.1 Inleiding 


Van 27 tot en met 30 april 1992 heeft in Krems (Oostenrijk) een congres plaatsgevonden 
onder de titel ‘Spring school on the didactics of computet algebra’. Het onderwerp was 
de didactische inpassing van computeralgebra, en dan met name Derive, in het voortgezet 
onderwijs. De organisatie was in handen van het Research Institute on Symbolic Computation 
(RISC) uit Linz. Voor de financién zorgde het Oostenrijkse ministerie van onderwijs terwijl 
Soft Warehouse Europe ondersteuning verleende. Vanuit het Freudenthal Instituut was ik bij 
deze conferentie aanwezig en hieronder volgt een impressie van het geheel. 


10.2 De voordrachten 


De deelnemers (zo’n 35 in totaal) waren afkomstig uit diverse landen: behalve natuurlijk 
Oostenrijkers en Duitsers waren er Spanjaarden, Britten, een Zweed, een Braziliaanse maar 
vreemd genoeg geen Fransen. Toch wel een beetje speciale deelnemers uit Hawai waren David 
Stoutemyer, de maker van Derive, en zijn vrouw. 

In principe hield elke deelnemer een voordracht. Veel daarvan hadden betrekking op ervarin- 
gen in een of meer klassen met de behandeling van een onderwerp met Derive. Al met al een 
veelheid aan uitleenlopende voorbeelden! Verschillende opvattingen over wiskundeonderwijs 
leiden natuurlijk ook tot een verschil in de wijze waarop Derive gebruikt wordt. Daarnaast 
zijn de curricula en de onderwijstradities in de diverse landen ook zeer verschillend. 

Toch waren er gemeenschappelijke elementen te signaleren. Derive wordt vaak gebruikt van- 
wege de combinatie van grafische en algebraische mogelijkheden. In een grafisch venster 
kan snel een visualisatie gemaakt worden, terwijl in een algebravenster het (symbolische) 
rekenwerk foutloos uitgevoerd wordt. Een voorbeeld van dergelijk gebruik uit een lezing van 
de heer Zochling over het gebruik van Derive bij het onderwijzen van de gaswetten in de 
natuurkundeles staat in figuur 1. 


Een tweede punt dat vaak terugkwam was de mogelijkheid om, nu Derive of een vergelijkbaar 
systeem het rekenwerk uitvoert, aandacht te besteden aan modelbouw, simulatie, patroon- 
herkenning, zelfontdekkend leren, probleemoplossen en meer van dergelijke mooie zaken.. 
Hierbij moet wel worden opgemerkt dat de meeste delenemers lesgeven op het hoogste nivo, 
zeg maar wiskunde B van het VWO, al was er ook een bijdrage over de stimulerende rol van 
Derive voor zwakkere leerlingen. 

Stoutemyer zelf verzorgde ook een workshop waarin 0.a. een geavanceerde manier van sub- 





1De vorige nieuwsbrief bevatte abusievelijk het artikel Derive in the Dutch Classroom van Paul Drijvers, 
dat bestemd was voor het conferentierapport van de 1992 Spring School on the Didactics of Computer Algebra. 
Het verslag van deze conferentie hebt u dus nog te goed en treft u hierbij aan. 
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stitueren met behulp van limieten aan de orde kwam. Een klein voorbeeld uit xijn workshop 
staat afgebeeld in figuur 2. Twee maal is de grafiek getekend van de geéxpandeerde vorm van 
(x — 1)°. Links is de grafiek grillig door afrondfouten. Rechts is de precisie verhoogd. Dat 
leidt natuurlijk tot een betere grafiek maar wat mij verraste is dat het tekenen met de hogere 
precisie ook nog sneller gaat dan bij de lagere nauwkeurigheid. Dat heeft te maken met het 
tekenalgorithme: als er door de afrondfouten veel variatie in de functiewaarden lijkt te zitten 
dan worden er meer punten uitgerekend wat tot vertraging leidt. Het lijkt paradoxaal: hogere 
precisie en tevens meer snelheid. 














FIGUUR 2. TWEEMAAL DEZELFDE GRAFIEK? 


10.3. Technologische ontwikkelingen 


Op de conferentie circuleerde een drietal palm-top computers van het type HP 95LX waarop 
Derive geinstalleerd is. Het is een apparaatje dat niet groter ia dan een calculator en dat 
opengeklapt kan worden. Derive versie 2.07 draait er zonder problemen op. Het zal duidelijk 
zijn dat er veel belangstelling voor dit speelgoed was. Nooit meer problemen met het re- 
serveren van een computerlokaal, geen tweetal studenten met maar één apparaat, en ook 
thuis dezelfde faciliteiten. ‘Personal technology’ noemen ze dat geloof ik. Naar mijn mening 
is een dergelijk apparaat een voorwaarde voor een werkelijke doorbraak van computeralgebra 
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in het voortgezet onderwijs. 

Verder is veel gesproken over wensen die men t.a.v. Derive heeft. Een verrrassing voor mij 
was de nieuwste versie van Derive, versie 2.10. Een belangrijke verbetering van de interface 
is het feit dat het splitsen van het scherm in een reken- en een tekenvenster geautomatiseerd 
is. Vanuit mijn ervaringen in de klas was een dergelijke verbetering zeer gewenst! Verder 
is 2.10 uitgebreider dan zijn voorgangers. Er is om maar eens wat te noemen een proce- 
dure DSOLVEI GEN die differentiaalvergelijkingen (voor zover mogelijk) oplost zonder dat 
de gebruiker eerst moet aangeven om welk type d.v. het gaat. Zo zijn er nog een aantal 
verbeteringen. 

Over een Derive-versie die met de muis te besturen is wordt momenteel hard nagedacht. Mo- 
gelijk verschijnt er een versie die onder Windows draait. 

Bij dit alles komt natuurlijk ook meteen de kwestie van licenties en prijzen aan bod. Het kan 
geen kwaad om in deze naar het buitenland te kijken. In Oostenrijk heeft de staat de levering 
van Derive aan alle gymnasia geregeld. Ook in Frankrijk staat iets dergelijks op stapel, zij 
het dat de Fransen als voorwaarde stelden dat Derive geen woord Engels bevat. In Engeland 
start binnenkort een project met een budget dat het mogelijk maakt om leerlingen van hun 
eigen HP 95LX te voorzien. Wordt het niet eens tijd dat ook in Nederland over een landelijke 
licentie wordt nagedacht? 


10.4 Een didactiek van computeralgebra? 


Of de titel van de conferentie waargemaakt is valt enigszins te betwijfelen. De didactiek 
van de computeralgebra lijkt vooralsnog te bestaan uit een stel mooie voorbeelden, leerzame 
ervaringen, hoopgevende intuities en vooral vragen die een nader onderzoek behoeven. Naar 
mijn mening is het duidelijk dat computeralgebra een belangrijk fenomeen in het wiskun- 
deonderwijs zal worden, maar zijn we nog lang niet toe aan een uitgekristalliseerde visie op 
hoe dat onderwijs er dan uit zou kunnen en moeten zien. 

In de slotdiscussie bleek eens te meer dat we allemaal hopen dat computeralgebra een bij- 
drage zal leveren aan een meer experimentele ‘bottom-up’ benadering van het leren van 
wiskunde. En dat er meer ruimte zal komen voor ‘hogere’ vaardigheden als modelbouw en 
mathematiseren en probleemoplossen, omdat het uitvoeren van algorithmen minder relevant 
wordt. Maar of dit allemaal te realiseren is, en of dit in het voordeel of in het nadeel van de 
zwakkere leerlingen is, en hoe een en ander didactisch vormgegeven moet worden, daarover 
bestaat nog veel onzekerheid. Ook de inhoud van het wiskundeonderwijs staat ter discussie, 
erkent iedereen. Maar als de vraag concreet wordt is de eensgezindheid ver te zoeken. Of 
leerlingen in de toekomst de regels voor het differentiéren nog moeten kennen? De meningen 
zijn verdeeld. 


10.5 En verder ... 


En verder was het vooral erg inspirerend en aangenaam om enkele dagen met ‘gelijkgestemde 
geesten’ ervaringen te delen, materiaal uit te wisselen, contacten te leggen enzovoorts. Verder 
hadden onze Oostenrijke gastheren voor een aantrekkelijk cultureel programma gezorgd, 
waardoor mijn Oostenrijks met het woord ‘heuriger’ is uitgebreid. Dat blijkt een soort wijn- 
café te zijn waar het goed toeven is. Wist u trouwens dat er een echte Derive Song bestaat, 
die door alle aanwezigen uit volle borst is meegezongen? Kortom, het was erg gezellig. 


Paul Drijvers 
Freudenthal Instituut 
Utrecht 
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11 | Computer Algebra Meetings 


Below is a list of upcoming computer algebra meetings like workshops, summer schools, 
conferences, etc. Note that from many meetings more information is available by anonymous 
ftp from canb.can.nl through the directory /ftp/pub/CA-meetings. 


1. Computeralgebra in het Technisch Onderwijs: Toverdoos of Black Box 
Noordelijke Hogeschool Leeuwarden, May 12, 1993 


Contact address: Ms. B. Vis, NHL Sector Techniek, Vondelstraat 9, 8913 HP Leeuwarden 


2. Helsinki Workshop on Symbolic and Numeric Computation 
University of Helsinki, Finland, May 24 — 28, 1993. 
Contact address: Anne Bjorklund, CSC-Tieteellinen laskenta Oy, PO Box 405, Fax: +358-0-4572302, Phone 
+358-0-4572158 

3. “Algorithm Week” in conjunction with “Groups and Geometries at O.S.U.” 


Ohio State University, Columbus, Ohio (USA), May 24 — 28, 1993. 


Contact address: Dr. A. Seress, akos@function.mps.ohio-state.edu 


4. Computational Commutative Algebra 
Cortona, Italy, May 31 - June 4, 1993. 
Contact address: L.Robbiano, Dipartimento di Matematica, Via L.B.Alberti 4, I-16132 Genova, email: 
robbiano@matgen.ge.cnr.it , fax: Italy + 10 + 3538769 
5. Konferenz der deutschsprachigen AXIOM-Benutzergruppe 
Karlsruhe, Germany, June 24 ~ 26, 1993. 
Contact address: Prof. Dr. G. Schneider, 0721-608-2479, Fax: 0721-32550, email: schneider@rz.uni-karlsruhe.de, 
schneider@dkauni2.bitnet; Dr. J. Grabmeier, 06221-594329, Fax: 06221-593500, email: grabm@dhdibm1.bitnet. 
6. Commutative Algebra and its Interaction with Computer Algebra and Com- 
binatorics 
Inst. f. Experimentelle Mathematik, Essen, Germany, June 28 —30, 1993. 
Contact address: Prof. Dr. J. Herzog, registration before end of May with Frau B. Ebinger, Inst. f. Exp. Math., 
Ellernstr. 29, 4300 Essen 12, email: mat425Q@vm.hrz.uni-essen.de, Tel. 0201-32064-52, Fax: -25. 
7. Maple Summer Workshop and Symposium 
Univ. of Michigan, Ann Arbor, June 28 — 30, 1993. 


Coordinator MSWS 93: Paola D’Alessandra, Waterloo Maple Software, 450 Phillip St., Waterloo, Ontario, 
Canada, N2L-5J2, Tel: (519) 747-2373, Fax: (519) 747-5284, email: msws@maplesoft.on.ca 


69 


70 


10. 


11. 


12. 


13. 


14. 


15. 
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International Conference on Number Theoretic and Algebraic Methods in 
Computer Science 


Moscow, Russia, June 29 — July 2, 1993. 


Conference chairman: Horst G. Zimmer, Saarbr” ucken, (zimmer@math.uni-sb.de) 
Contact address: Igor Shparlinski, Dept. no. 4, ICSTI, Kuusinenastr., 21 B, Moscow, 125252, Russia FAX 
(095) 943 00 89, Email: shparplb.icsti.su, 


ISSAC’93 
Kiew, Ukraine, July 6 — 8, 1993. 


Contact address: Alexander A.Letichevsky, EECS Dept., Comput. Sci. Div., 571 Evans Hall, Univ. of California, 
USA-Berkeley, CA 94720, email: fateman@cs.berkeley.edu, +1-510-642-1879 Programme committee: Gaston 
Gonnet, ETHZ, Zentrum, CH-8092 Ziirich, email: gonnet@inf.ethz.ch. 


NSF Workshop for Factuly with Symbolic Experience 
Rose-Hulman Institute of Technology, USA, July 13 - 17, 1993. 


Contact address: Rose-Hulman Institute of Technology, 5500 Wabash Avenue, Terre Haute, IN 47803, fax: 
(812) 877-8895, email: NSF.Workshop@Rose-Hulman.edu, staff: Dr. Mark A Yoder, Electrical Engineering, 
phone (812) 877-8291, email: Mark.A.Yoder@Rose-Hulman.edu Dr. Robert Lopez, Dept. of Mathematics, (812) 
877-8396, email R.Lopez@Rose-Hulman.edu 


IMACS Symposium on Symbolic Computation ’93 
Lille, France, July 14 -17, 1993. 


Contact address: G. Jacob, France (chairman), Symposium IMACS SC 93, Bat. M3-Informatique Université 
Lille I, F-59655 Villeneuve d’Ascq Cedex France, email: jacob@lifi.fr, Stan Steinberg, U.S.A. (co-chairman) 
email: stanlyQmath.unm.edu. 


MTNS 93: Computer Algebra in Systems Theory and Control 
Regensburg, Germany, August 2 — 6, 1993. 


Contact addres: M. Hazewinkel, CWI, Kruislaan 413, 1098 SJ Amsterdam, The Netherlands, email: mich@cwi.nl 


Summerschool “Computational Group Theory” 

Galway, Ireland, August 1 -14, 1993. 

Contact address: Dr. E. Robertson, Groups 1993, Mathematical Institute, Univ. of St. Andrews, KY16 9SS, 
Scotland 

Cayley/MAGMA Conference on Computational Algebra 

London, England, August 23 - 27 1993. 


Contact address: R.L. Siverstone, School of Mathematical Sciences, Queeen Mary and Westfield College, Mile 
End Rd., London El 4NS, England, email: R.L.Siverstone@qmw.ac.uk. 


Design and Implementation of Symbolic Computation Systems DISCO’93 
Gmunden, Austria, September 15 — 17, 1993. 


Contact address: Dr. Franz Lichtenberger, RISC, Johann Kepler Universitat, A-4040 Linz, Austria, email: 
lichtenbGrisc.uni-linz.ac.at Programme committee: Prof. Alfonso Miola, DIS Universita’ di Roma “La Sapienza” 


Via Salaria 113, 00198 Roma, Italy, email: disco93@discol.ing.uniromal.it 
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16. Technology in Mathematical Teaching 
Birmingham, England, September 17 — 20, 1993. 
Contact address: Pam Bishop, CTICMS, Fac. of Education, The University of Birmingham, Edgbaston, UK- 
Birmingham, B15 2TT. Email: tmt93@bham.ac.uk 

17. Teaching Mathematics with Derive 
Krems, Austria, September 25 — 28, 1993. 
Contact address: RISC-Linz/Krems’93, Attn. Bernhard Kutzler, University of Linz, A-4040 Linz, Austria, fax: 
+43-7236-3231-30 

18. Workshop on Software for High Energy and Nuclear Physics 
Oberammergau, Germany 


Contact address: Verena Offermann, Physics Department, University of Wuppertal, PO Box 100127, D-5600 
Wuppertal 1, Germany, Tel: +49-202-439-2733, Fax: +49-202-439-2811, 
E-mail offerman@physik.uni-wuppertal.dbp.de 


19. Symposium on Quantifier Elimination and Cylindrical Algebraic Decompo- 
sition 
RISC Linz, October 6 — 8, 1993. 


Contact address: Prof. Jeremy R. Johnson, Department of Mathematics and Computer Science, Drexel Univer- 
sity, Philadelphia, PA 19104 USA, Tel.: 001-215-895-2893 email: jjohnson@mcs.drexel.edu. 


More information: Prof. B. F. Caviness, Department of Computer & Information Sciences, 103 Smith Hall, 
University of Delaware, Newark, DE 19716 USA, Tel.: 001-302-831-8234 email: cavinessQudel.edu. 


20. MEGA 94: Effective Methods in Algebraic Geometry 
Santander, Spain, April 5 — 9, 1994. 


Contact address: Prof. Tomas Recio, Prof. L. Gonzdlez—Vega, Deadline for contributions is 16.09.1993 in TeX 
in electronic form to mega94@ccucvx.unican.es or 16 copies to Departamento de Matematicas, Facultad de 
Ciencias, Universidad de Cantabria, Santander 39071, Spanien, Tel.: 0034-42-201433, Fax: 0034-42-201402, 


email: recio@ccucvx.unican.es. 


CAN account 


subscription form 


ref: accl1993 


CAN Software and Computers 


The CAN Expertise Centre has computing facilities that are dedicated to computer algebra. These 
facilities are accessible via SURFnet, Internet and Datanet. A series of computer algebra systems are 
installed on the machines: general purpose systems like Aljabr (Macsyma), FORM, Maple, Mathematica, 
Reduce; and also special purpose systems like DStool, GAP, GB, LiE and Macaulay. CAN tries to 
get as many computer algebra systems and packages on these machines as possible. At this moment 
the accessible hardware consists of a multiprocessor Silicon Graphics Crimson and two Silicon Graphics 
Indigo’s with R4000 processors. 


Conditions 


1. The private disk space for the user is limited to 5 MB. 

2. The account can be used up to December 31, 1993. 

3. The use of CAN computing facilities is strictly personal, the password must be kept secret. Abuse will 
lead to a denial of further access. 

4. The use of the CAN computing facilities is regulated by the system operators. They are allowed to 
stop processes, or to refuse access, if necessary. Before very large computations on the CAN computing 
facilities, the systemoperator should be notified. 

5. This subscription is governed by the Dutch law. 


Subscription Form 


I hereby ask for an account in 1993 on the CAN computers for the price of 350 DAl. exclusive VAT 
(app. 210 US$). I confirm to have read and to agree with the conditions above. 























name: 
company /institute: 
address: 

city: postal code: 

country: e-mail address: 

tel: fax: 








desired login name: 





method of payment: 


O send bill to above address 


O check is enclosed; 





O charge VISA cardno___—————S—_C*7"= expiration date 


date: ___—=————CC:C Signature: 
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THE COMPLETE SOLUTION 
FOR TECHNICAL COMPUTING 














Visualization 











When it was released in 
| 1988, Mathematica was hailed 
as a software system that 
“fundamentally alters the 
mechanics of mathematics.” 
Four and a half years later, 
with several hundred thou- 
sand scientists, engineers, 
students, financial ana- 








Macintosh 














Algebra and Caiculus 
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lysts, medical researchers, 
and others using the 
system, Mathematica is 
firmly established as the 
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Simplify [D[%, x]] 





| | log{x_ y_] := log{x] + logfy] 
© Numbers F 


log[a b c*2 4d] 
Ini] | | ower 


N[Sqrt [Pi], 250] 


2453850905516! | log[x_ * n_] := n log[x] 
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738544665416226| | Logfa b c*2 £) 


236152865724422| | owe 
787068462037698|,] log[a] + log{b] + 2 


log{a] + log{b] + 1 


a ||| standard environment 
| for technical computing. 








x x 5Symbolic Programming d Starting Out 


When you first get start- 
ed, you can use Mathe- 
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cate} + logtai | matica just like a cal- 


culator. But it is a calcu- 
lator that lets you do cal- 
Log Let +: Lemke? | culations as large as you 
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© Graphics 
Info} = 
DensityPlot [Sin[Re[Gamma[x+I y]]], 
{x, 3, 7}, fy, 0, 5), 
PlotPoints->200, Mesh->False] 























Numerical computation with unlimited precision integer, real, and complex numbers. Exact symbolic computation. Wide 
range of mathematical functions, including hypergeometric special functions, number theory functions, statistics. Nu- 
merical and symbolic matrix operations. Numerical and symbolic equation solving, root finding, integration, differential 
equations, optimization. 2D and 3D color graphics, animation, and sound generation. Graphics language for representa- 
tion of arbitrary structures. Output in standard PostScript. Data importing and exporting. External program linkages. Full 
built-in symbolic programming language, with procedural, functional, and rule-based capabilities. Notebook front end 
with hierarchical documents, style sheets, importing and exporting of graphics, sound, animations, and typesetting. Net- 
| work operation supported. Journal, newsletters, more than 30 books on Mathematica available. Add-on packages; free 
MathSource™ electronic resource. Versions: Microsoft Windows * Macintosh ¢ MS-DOS ¢ NEC PC # SPARC © DEC 
RISC, VAX # HP ¢ IBMRISC # NeXT © SGI © Convex ¢ and others. Student versions available. Now shipping Version 2.2. 


— want, and that can han- 
dle the full range of mathematical 
functions, matrices, algebraic ma- 
nipulation, and calculus as well as 
creating publication-quality 2D and 
3D graphics. 


Flexible Programming 
The easy-to-use language that is 
built into Mathematica requires no 








learning for simple calculations, 
but has the power to let you create 
even the most sophisticated pro- 
grams. And with the Notebook con- 
cept pioneered by Mathematica, 
you can create interactive docu- 
ments that mix text, graphics, ani- 
mations, and sounds with live 
formulas and programs. 


A Range of Applications 
Whether you design airplanes, solve 
quantum field theories, analyze 
medical experiments, or just want to 
do your algebra homework, Mathe- 
matica provides a unique environ- 
ment for solving your problems. 
And you can run your same Mathe- 
matica programs without change on 
everything from Macintosh and PCs 
(with or without Windows) to Unix 
workstations and supercomputers. 


Wolfram Research 


Wolfram Research, Inc. 
100 Trade Center Drive, Champaign, IL 61820-7237, USA 
+1-217-398-0700; fax: +1-217-398-0747; email: info@wri.com 


Wolfram Research Europe Ltd. 

Evenlode Court, Main Road, Long Hanborough, Oxon OX8 2LA, UK 
+44-{0)993-883400; fax: +44-(0)993-883800 

email: info-euro@wri.com 


Representatives in over 30 countries; contact main office 





©1993 Wolfram Research, Inc Mathematica is a registered trademark and MathSource is a trademark of Wolfram Research, inc. Mathematica is not associated with Mathematica Inc., Mathematica Policy Research, Inc. or MathTech, Inc. All other product names mentioned are trademarks of their producers. 





